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important  in  plasmas  of  appreciable  fractional  ionization. 


These  effects  can  be  adequately  described  by  the  Fokker- 
Planck  treatment.  Unfortunately,  the  resulting  expression 
greatly  complicates  the  analysis  and,  for  this  reason,  is 
often  ignored.  In  this  investigation,  the  Fokker-Planck 
expression  is  replaced  by  a simple  relaxation  term. 
Although  this  approach  introduces  considerable  error,  it 
does  serve  to  indicate  several  important  features  which 
must  be  included  to  represent  the  physical  situation 
adequately.  Hopefully,  this  information  will  be  useful 
in  any  further  attempts  to  simplify  the  Fokker-Planck 
treatment . 

I am  grateful  to  Maj.  P.  E.  Nielsen  and  Capt. 
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Abstract 

' 

The  electron-electron  collision  term  of  the  steady- 
state  Boltzmann  equation  is  replaced  by  a simple  relaxation 
term  that  effectively  linearizes  the  equation.  Additional 
assumptions  are  made  to  simplify  the  equation  further. 

Ionization  is  ignored  and  a two-level  atom  is  assumed.  The 
effects  of  a DC  electric  field  are  included.  Using  the 
relaxation  term  to  account  for  electron-electron  collisions, 
an  approximate  analytic  solution  is  derived.  Temporal  and 
steady-state  characteristics  of  the  relaxation  term  are 
compared  to  those  of  the  standard  Fokker-Planck  term.  The 
relaxation  term  is  judged  invalid  for  energies  greater  than 

■ 

the  excitation  threshold. 


I Introduction 


There  exists  a large  class  of  problems  in  plasma 
physics  where  important  physical  quantities  are  sensitive 
to  the  detailed  form  of  the  electron  distribution  function. 

To  predict  these  quantities  accurately,  the  time-independent 
Boltzmann  equation  must  be  solved. 

This  equation  is  very  difficult  to  solve  in 
general.  Some  simplifying  assumptions  must  be  made,  even 
to  solve  the  equation  numerically.  In  addition,  terms 
are  often  neglected  on  the  basis  of  complexity.  The 
effects  of  electron-electron  collisions  are  particularly 
difficult  to  handle  because  of  the  non-linear  terms  intro- 
duced , 

For  high  electron  number  densities,  electron-electron 
collisions  can  be  very  important.  These  interactions  tend 
to  drive  the  distribution  toward  a Maxwellian  and  can 
appreciably  alter  the  high  energy  part  of  the  distribution. 
Therefore,  there  exists  a need  for  a simple,  but  accurate, 
term  to  describe  the  effects  of  electron-electron  collisions. 
In  an  effort  to  fulfill  this  need,  a simple  Krook-type 
relaxation  term  is  developed  and  its  characteristics  are 
investigated. 

The  theory  leading  to  the  Boltzmann  equation  and 
the  Fokker-Planck  collision  term  are  reviewed  in  Section 
II.  With  this  background,  a simple  relaxation  term  is 
developed  and  an  approximate  solution  to  the  Boltzmann 
equation  is  derived  in  Section  III.  The  results  of  the 


1 


7 3raRBK3CBaHBR?27*  t*" 


relaxation  term  and  the  Fokker-Planck  term  are  analyzed 
and  compared  in  Section  IV.  Finally,  the  conclusions 
reached  and  recommendations  for  future  work  are  presented 
in  Section  V. 
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II  Theory 


In  this  section,  the  problem  is  carefully  defined 
and  the  governing  equation  is  derived.  The  assumptions 
necessary  to  arrive  at  this  equation  are  stated  as  they 
occur  in  the  derivation.  The  macroscopic  effects  and 
non-Coulomb  interactions  are  incorporated  using  the 
classical  Boltzmann  formalism.  The  Coulomb  interactions 
are  included  using  the  Fokker-Planck  approximation. 

Finally,  the  assumptions  are  justified  using  order  of 
magnitude  estimates  where  possible. 

The  Boltzmann  Equation 

The  plasma  considered  here  is  assumed  to  consist 
of  electrons,  singly-charged  ions,  and  neutral  particles. 
Internal  excitation  levels  are  ignored  for  the  ions 
and  only  one  level  (above  the  ground  state)  is  allowed 
for  the  neutrals.  Ionization  is  ignored  for  both  ions 
and  neutrals.  Electron  loss  mechanisms,  such  as  attach- 
ment and  recombinations,  are  also  ignored.  Therefore, 
the  number  densities  of  the  various  particles  can  be 
treated  as  parameters  in  what  follows. 

It  is  assumed  that  the  equations  of  kinetic 
theory  are  applicable.  This  is  acceptable  provided  (Ref  13) 


Kkt*) 


•/i 


»/* 


» * 
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where 


m*  = mass  of  species  <* 

T*  = kinetic  temperature  of  species  <* 

= number  density  of  species  <* 

K = Boltzmann's  constant 
= Planck's  constant/27T 

In  accordance  with  kinetic  theory,  a velocity  distribution 
function  (r ,v,t)  is  defined  for  each  species  such  that 
the  average  number  density  of  particles  of  type  <x  with 
velocities  between  v and  v + dv,  located  in  the  spatial 
volume  element  dr  centered  at  r,  at  time  t is  given  by 
F^ (r,v, t)drdv>  the  average  being  carried  out  over  a macro- 
scopically  infintesimal  time  interval.  Ground  state  neutrals 
and  excited  neutrals  are  to  be  regarded  as  separate  species. 
The  distribution  functions  are  determined  by  the 

Boltzmann  equation.  This  equation  may  be  written  symbolically 
as 


Pfftr  / i O — ^ FvCr, 

cD{F<(E,v/t)}  - ,co(j 


(1) 


where 


c)*t 


(2) 


and 

Fa  = external  macroscopic  force  on  species 
Vv=  gradient  operator  in  velocity  space 

= gradient  operator  in  configuration  space 


4 


represents  the  change  in  F^  resulting 


J R^r.a.Ol 

» c»!| 

from  the  interaction  between  particles.  Macroscopic 

effects,  such  as  external  forces,  are  usually  included 

in  JD" . Microscopic  effects,  such  as  excitation  by 

electron  impact,  are  usually  included  in  ^€.Ly>^-|  . 

o-t  »C*|| 

In  general,  the  change  in  F^  resulting  from  interactions 
with  other  species  depends  on  the  distribution  functions 
of  these  species.  Therefore,  this  term  couples  the 
equation  for  one  species  to  all  the  other  equations. 

Fortunately,  for  many  problems  a knowledge  of  the 
electron  distribution  is  sufficient  to  describe  the 
important  features  of  the  system  adequately.  To  simplify 
the  problem,  this  assumption  is  made  here.  Furthermore, 
the  distributions  of  the  other  species  are  assumed  Maxwellian 

With  these  assumptions,  a complete  statistical 
solution  requires  that  the  electron  distribution  function 
F(r,v,t)  be  determined  using  Eq  (1).  Note  that  F(r,v,t), 
with  no  subscript,  will  be  used  exclusively  for  the  electron 
distribution  function. 

It  is  assumed  that  F(r,v,t)  is  spatially  homogeneous 
and  the  only  external  forces  result  from  a constant  electric 
field  E . Since  F is  now  independent  of  r,  it  will  be 
denoted  by  F(v,t).  The  second  term  of  Eq  (2)  vanishes 
in  this  case  and  Eq  (1)  becomes 

- li . y f(v  t)= 

Vvr^-'  ; a*  lCol, 


(3) 


where 


m = electron  mass 
-e  = electron  charge 
ne  = total  number  density  of  electrons 


The  normalization  is  such  that 


■V  = j F(v, 


O d V 


In  general,  nQ  should  be  expressed  as  a function  of  the 
time,  but  the  neglect  of  electron  sources  and  sinks  requires 
that  ne  be  constant. 

The  right  side  of  Eq  (3)  involves  the  interactions 
of  the  electrons  with  themselves  and  the  other  plasma 
species.  It  is  assumed  that  these  terms  impart  little  net 
motion  to  the  electrons.  In  fact,  these  interactions  should 
tend  to  produce  an  isotropic  velocity  distribution.  The 
second  term  in  Eq  (3)  accounts  for  the  external  field  effects. 
This  term  does  produce  anisotropic  effects.  If  the  field  is 
not  too  large,  the  electrons  will  be  accelerated  very  little 
between  the  randomizing  collisions.  The  principal  effect  of 
the  field  is  a heating  of  the  electrons;  the  anisotropic 
drift  in  the  direction  of  the  field  is  a relatively  minor 
effect. 

The  above  arguments  motivate  an  expansion  of  F(v,t) 
in  spherical  harmonics.  For  small  anisotropic  effects  a 
first  order  expansion  is  adequate  and  F(v,t)  is  given  by 


F0(v,-t)  + F,(v,t)  cos  0 
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where  © is  the  angle  between  v and  the  applied  field  E 
If  Eq  (5)  is  substituted  into  the  left  side  of  Eq  (3 ) » the 
result  is 


+ cos cose 

o£  rn  £t 

eE  . 3F,(v,0  e£  $\nz&  r- 1 .j.\  — I ( 

" l?rtos  0 ~rt ~ - v—  Kw.t)  - at  lc«M  ' } 

Equation  (6)  is  now  integrated  over  all  solid  angles  d-H-  , 


where 


d-Ti-—  Sirt  © d 


and  <$)  is  the  azimuthal  angle  about  the  polar  axis.  The 
resulting  equation  is 


^ - Is  !>  h ^ »>^J  = SF  aiL  (7) 


If  Eq  (6)  is  multiplied  by  cos  & and  again  integrated  over 
all  solid  angles,  Eq  (8)  is  obtained* 


^F  (vt)  _ e.  E-  ^ F0(v,t) 
<it  nn  & v 


3 F(y,0 
Cos  9 at 


The  collisional  integrals  in  Eqs  (7)  and  (8)  remain 
to  be  evaluated.  The  analysis  can  be  simplified  by  consi- 
dering each  process  separately.  Since  electron  sources  and 
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sinks  have  been  ignored,  the  collision  term  takes  the  formi 


aFtv.O  _ 

Colt 


* lift 


^FCv.U 

b-t 


tn 

Uelas^ic. 


b FCv.t) 


ent 

elastic 


^F(y,Q 

St 


i**«las4-i«. 


3F(y,0 

— 5T~ 


ei 


<5F(v,0 

St 


eC 


(9) 


The  first  term  on  the  right  side  of  Eq  (9)  represents 
the  effect  of  elastic  collisions  with  neutral  particles  in 
the  ground  state.  If  the  neutral  particles  are  assumed 
infinitely  massive,  these  collisions  will  have  no  effect  on 
an  isotropic  distribution.  Therefore, 


•H 


= o 


For  an  anisotropic  distribution,  these  collisions  will  tend 
to  reduce  the  degree  of  asymmetry.  Using  the  standard  form 
of  the  Boltzmann  collision  integral,  it  can  be  shown  that 
(Ref  23*82) 


3F,(v,U 


en  — 
classic. 


where  V*,  is  known  as  the  momentum  transfer  collision 
frequency  for  electron-neutral  collisions.  The  collision 
frequency  may  be  expressed  in  terms  of  the  momentum  transfer 
cross  section  3^(v)  * 


■‘’i,  = %v<rm(v) 

where 

no  » total  number  density  of  neutral  particles 
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In  general,  is  a function  of  the  electron  velocity. 

However,  it  is  assumed  here  that 


= 


where  R is  a constant  momentum  transfer  rate.  Therefore 
m 


and  (from  Appendix  D),  Rm  = 5*928  x 10“®  cm^-sec”*.  With 
this  choice  for  R^,  the  collision  frequency  used  here  is 
approximately  the  same  as  that  for  electrons  and  ground  state 
Argon  atoms. 

The  second  term  in  Eq  (9)  represents  the  effect  of 
inelastic  collisions  (excitation)  with  neutral  particles  in 
the  ground  state.  These  interactions  convert  high-energy 
electrons  to  low-energy  electrons.  Inelastic  collisions 
have  a large  effect  on  the  bulk  of  the  distribution,  but  a 
minor  effect  on  anisotropities  (Ref  12).  Therefore 


J)  F,  (v,t) 


en 


— O 


The  effect  on  FQ(v,t)  can  also  be  calculated  using  the 
Boltzmann  collision  integral  formalism,  but  the  result  is 
most  conveniently  expressed  in  terms  of  energy.  For  now, 
this  effect  will  be  denoted  by 


a Rs(v,-Q 

at 


en 

»*el*vV‘c. 
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The  third  and  fourth  terms  in  Eq  (9)  represent  the 


above  effects  for  excited  state  neutrals  as  opposed  to 
ground  state  neutrals.  These  interactions,  known  as 
superelastic  collisions,  convert  low-energy  electrons  to 
high-energy  electrons.  These  terms  are  neglected  by  assuming 
negligible  excited  state  populations.  This  assumption  is 
discussed  further  at  the  end  of  this  section. 

The  fifth  term  represents  the  effects  of  the  Coulomb 
interactions  between  ions  and  electrons,  and  the  sixth 
between  electrons  and  electrons.  Since  the  Boltzmann 
collision  integral  is  based  on  short-range,  binary  collisions* 
this  formulation  is  not  well  suited  to  the  long-range  Coulomb 
interactions.  The  Fokker-Planck  equation,  presented  in  the 
next  section,  is  appropriate  in  this  case.  The  Fokker-Planck 
equation  will  be  used  only  for  the  electron-electron  collision 
term.  The  results  for  the  electron-ion  collision  term  are 
relatively  simple  and  will  only  be  quoted.  If  the  ions  are 
assumed  infinitely  massive,  the  argument  presented  for 
elastic  collisions  with  neutrals  yields 


^ F»(v,Q 

at  et  - 


O 


It  can  be  shown  (Ref  23«82) 


2>F.UO 

a-t 


~Ve.  i R C'/jO 
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where  is  an  effective  electron-ion  collision  frequency. 

It  is  shown  at  the  end  of  this  section  that  for 

the  fractional  ionizations  considered  here.  Therefore, 
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electron-ion  collisions  will  be  ignored  entirely.  In 
addition  will  be  ignored.  Electron-electron 

ot.  let 

collisions  will  tend  to  reduce  the  anisotropity  of  the 
distribution.  However,  for  the  fractional  ionizations 
considered  here,  the  elastic  collisions  with  neutrals  are 
the  dominant  effect  (Ref  21i284).  It  remains  to  calculate 
j , The  next  section  is  devoted  exclusively  to 
determining  this  term. 

The  collision  term  may  now  be  expressed  as 

aF(v,Q]  JF„(v,tl|  ,,, 

» U —St—  Li + c“ s >*L  • Ll  (lt 


where 


f (v  O 

AFo^l  r ar.tv.oi  _ . 

at  L,  Ll\+,c  + -TE — 


Using  Eqs  (10)  through  (12)  in  Eqs  (7)  and  (8)  results  in  a 
set  of  coupled  equations  independent  of  Q and  <P  , The 
steady-state  results  are 


fkili  ^ 3 R.  (v\ 

^t  iite-Ustic- 


_<JE  dF„U)  _ 
m J v 


- F,L\/) 


Equation  (14)  is  now  substituted  into  Eq  (13 ) • The 
resulting  time-independent  equation  involves  only  F0(v), 
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and  if  Vm  is  independent  of  V the  desired  expression  is 


e*  E* 
(p**1  vm 


*F.<v)| 

vjT  ] ft*'  . -V 
'loddsii* 


ee 


(15) 


Equation  (15)  may  be  transformed  from  velocity  space  to 
energy  space  in  a straightforward  manner.  Let 

and 


The  number  density  of  electrons  with  energies  between  £ 
and  £+<*£  is  n(£)d£  . Equation  (15) » after  the  appropriate 
transformations,  takes  the  following  formi 


dnCc)  7 _ ^ 

TT  i - . 


flu 

«€las4ic. 


dr>(£) 

^t.  ee 


(16) 


where 


(17) 


The  second  term  in  Eq  (16)  may  be  deduced  from  physical 
arguments.  A more  rigorous  treatment  is  given  by  Nielsen 
(Ref  18).  Let  x denote  the  excitation  energy  of  the 
neutral  particles,  and  R(£)  the  rate  at  which  excitations 
occur  for  electrons  with  energy  £ . The  excitation  rate  is 
related  to  the  total  excitation  cross  section  <TX  by 


(18) 
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Inelastic  collisions  can  affect  the  electron  distribution 
function  in  two  ways.  Firstly,  electrons  with  energy  £, 
can  produce  excitations  at  the  rate  R(£)  and  lose  energy 
x in  the  process.  This  represents  a net  loss  of  electrons 
with  energy  £ . Clearly,  this  loss  is  proportional  to  the 
number  density  of  electrons  with  energy  £ and  the  number 
density  of  neutral  particles  in  the  ground  state.  Secondly, 
electrons  with  energy  £+x  can  produce  excitations  at  the 
rate  R(£+x)  and  lose  energy  x in  the  process.  This 
represents  a net  gain  of  electrons  with  energy  £.  . This 
gain  is  proportional  to  the  number  density  of  electrons  with 
energy  £ + x and  the  number  density  of  unexcited  neutral 
particles.  Therefore,  the  second  term  in  Eq  (16)  may  be 


expressed  as 


d n(£^ 
at 


It n - n0R(£+xWf+x)-  noR(£)h(0 


As  previously  mentioned,  superelastic  collisions  are  ignored 
in  Eq  (19).  The  form  of  R(£)  depends  on  the  form  of  <rx(£>  . 
The  determination  of  0^(0  is  of  great  practical  interest, 
but  will  not  be  considered  here.  It  is  assumed  here  that 

<T*(0  = R(£)  (20) 

where 


R(e)  - 


£ < X 


tax 
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From  Appendix  D 


R0  * 1.703  x 10“8  cm^  - sec  (102) 

With  this  choice  for  Rq,  the  excitation  rate  used 
here  is  approximately  equal  to  that  for  excitation  from  the 
ground  state  of  Argon.  Equations  (16)  and  (19)  may  be 
combined  to  yield  the  following  form  of  the  steady-state 
Boltzmann  equation 1 

IF  ^ «U.[_  zT  ” ^ = R CE+v')  n„R(£)n(£)+-  ^£*)ee  (22) 

where  i is  given  by  Eq  (1?)  and  R(t)  is  given  by  Eq  (21). 
The  explicit  form  of  |ee  is  deduced  below. 


The  Fokker-Planck  Elec tr on-Elec tron  Collision  Term 

To  derive  the  expression  for  the  electron-electron 
collision  term,  it  is  necessary  to  digress  temporarily  and 
consider  the  electron  velocity  distribution  F(v,t).  From 
Appendix  A,  the  effect  of  electron-electron  collisions  may 
be  described  using  the  Fokker-Planck  equation  in  the 
following  form 1 


2-t 


+ j.  y b}: 

2-  A- 


F(y 


where 


<AVi>  — ~ 4>(  J/jA  V.)  d ( a*\ 

<AviAvJ>=  Av/j  nv.AirtdUy) 


(58) 


(56) 

(57) 
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and  is  the  probability  that  an  electron 

changes  its  velocity  from  y-Ay  to  y in  a time  A-t  as  a 
result  of  electron-electron  collisions.  For  an  infintesimal 
time  interval  at,  may  be  expressed  in  terms  of 

the  Coulomb  differential  scattering  cross  section  <T~  as 
(Ref  22) 

(23) 

where 


A U 

M. 

c : 


v - v' 

I VI  - v'{ 


( e* 

r 1 1 

\2  + ) 

Lsmo/oJ 

and  0 is  the  scattering  angle  in  tho  center-of-mass 
coordinate  system.  Using  Eq  (23),  Eqs  (56)  and  (57)  can 
be  rewritten  as 

<av->  = JJa  vc  F(  v' t)  <r(x±,sz)  dy'd-fL  (24) 

^Av^AVj-)-  JJav/jAV;  F(  y'-t)a-(^Jll')  d v'dil-  (25) 

The  above  integrals  are  easier  to  calculate  in  a local 
coordinate  system  oriented  with  respect  to  the  relative 
velocity  w between  the  scattering  electrons.  It  is  worth 
noting  that  w is  not  the  velocity  change  Ay  of  a parti- 
cular electron  upon  scattering.  However,  for  elastic 
scattering  between  particles  of  equal  mass  Ay  = £ uy . 
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A fundamental  problem  arises  when  performing  the  above 
integrations!  the  integrals  diverge  logarithmically  at 
small  scattering  angles.  This  is  due  to  the  long-range 
nature  of  the  Coulomb  force.  It  can  be  shown  that  the 
Debye  length  is  the  maximum  distance  over  which  mutual 
electron  interactions  are  effective  (Ref  5)«  This  provides 
a natural  limit  in  performing  the  integrations;  scattering 
between  electrons  with  impact  parameters  greater  than  the 
Debye  length  is  ignored. 

If  the  integrals  are  evaluated  in  the  local  coordi- 
nate system  and  transformed  back  into  the  fixed  coordinate 
system,  the  result  is  (Ref  22) 

<Av->=r!£  <26) 

«W;av/,.>  =p£i_  (27) 

2>VC  ^ v • 

J 

where 


h r 2 

J'  FC 

(28) 

*• = T 

(29) 

r - jt»  a 

Vv\ 

4 

(30) 

/ 1 \'lz. 

(31) 

2 e3  \ rr  ne  J 

In  Eq  (31) » KT  is  proportional  to  the  average  electron 
energyi  for  no n-Waxwe Ilian  distributions,  T can  be 
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regarded  as  an  effective  electron  temperature.  The  quantity 
A is  actually  the  ratio  of  the  Debye  length  to  the  classical 
distance  of  closest  approach  for  two  electrons  with  relative 
velocity  xc  . Since  A appears  as  the  argument  of  a loga- 
rithm, it  is  customary  to  neglect  this  velocity  dependence 
and  use  an  average  relative  velocity  = J3  KT/m 

In  principle,  Eqs  (26)  through  (29)  can  be  substituted 

into  Eq  (58)  to  yield  an  expression  for  1 . The 

o"t.  1^^ 

above  equations,  despite  their  simple  form,  present  consider- 
able computational  difficulties  for  an  arbitrary  distribution 
F(v,t).  However,  for  the  isotropic  distributions  considered 
here,  these  results  may  be  further  simplified.  The  result 
for  an  isotropic  distribution  FQ(v,t)  is  (Ref  19) 


- rfi  o..o a.  Wu*] 

+ Fl(x.O  + jj-  [ JxF(«,t)dx 

o 

-/*F(x,^o-$-ro+&Oci*]]  <32> 


Equation  (32)  may  be  transformed  from  velocity  space  to 
energy  space  using  the  same  procedure  as  for  Eq  (15).  The 
resulting  equation  is  conveniently  expressed  in  the  following 
"flux-divergent"  form  (Ref  19) i 

^nC£«t)  | _ "b  ) 

"It  lee 


where 

T ~ ^[p(rr~  "ft) - ^ (3*0 
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and 

P = Z Z'',X  fxn(x,t)Jx  + 2£  jV'SKOdx 

« £ 

Q = 3 ^‘,/l  f »(*,*)  dx 

6 

« = frreV^)'^A 

Equation  (33)  is  the  final  form  of  the  electron-electron 
collision  term  to  be  used  in  Eq  (22).  Since  Eq  (22)  is  a 
steady-state  equation,  the  time  dependence  in  Eqs  (33) 
through  (3 6)  can  be  suppressed. 

A Discussion  of  the  Assumptions 

The  assumptions  made  in  deriving  Eqs  (22)  and  (33) 
deserve  closer  examination.  These  assumptions  can  be 
categorized  as  follows* 

(a)  Idealizations  have  been  used  to  model  the 
problem  in  a simple  fashion. 

(b)  Effects  which  are  small  compared  with  the 
dominant  effects  have  been  ignored. 

Assumptions  of  the  first  type  cannot  be  expected  to 
be  valid  in  general.  The  problems  considered  must  be  re- 
stricted to  those  with  essentially  the  same  physical  charac- 
teristics as  the  model.  Assumptions  of  the  second  type  can 
be  rationalized  using  physical  arguments  or  justified 
quantitatively  if  the  magnitudes  of  the  various  effects  can 
be  estimated. 

The  systems  considered  here  are  restricted  to  DC  gas 
discharges  such  as  those  commonly  used  to  pump  lasers.  The 


(35) 

(36) 
(3?) 
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field  is  assumed  uniform  and  there  is  no  mass  motion  of  the 

gas.  For  these  conditions,  the  electron  distribution 

function  can  be  considered  spatially  homogeneous.  Only  one 

neutral  ground-state  species,  with  a number  density  of 

lO^^cm”^,  is  considered.  This  is  closely  approximated  by 

a number  of  actual  laser  systems,  where  one  species  accounts 

for  over  90  per  cent  of  the  total  gas  mixture,  and  the  number 

density  chosen  is  also  representative  of  many  actual  laser 

systems  (Ref  15)*  The  applied  electric  field  E may  now  be 

expressed  in  terms  of  E/ne  , where  E is  measured  in 

volts/cm  and  nQ  is  the  number  density  of  the  neutral 

ground-state  species.  The  range  considered  here  is 

— 1 ft  —1  O 

restricted  to  10“  to  10“  volts-cm“.  This  range  is  well 
within  the  stable  discharge  region  for  most  systems  (Ref  15). 
Detailed  numerical  calculations  for  several  actual  laser 
systems  indicate  that  a negligible  amount  of  energy  goes  into 
ionization  for  the  E/r>0  range  considered  here  (Ref  20). 
Electron  attachment  and  recombination  are  the  most  significant 
electron  loss  mechanisms.  The  cross  sections  for  these 
processes  vary  depending  on  the  species  involved  (Ref  29«4l). 
If  attachment  or  recombination  is  significant,  it  is  neces- 
sary to  include  the  effects  of  ionization  to  insure  an 
appreciable  electron  density.  The  cross  sections  for 
multiple-ionization  are  usually  zero  until  very  high  energies 
and  then  orders  of  magnitude  lower  than  that  for  single- 
ionization. This  effect  may  be  safely  ignored  in  almost  all 
calculations. 

From  the  preceding  discussion,  it  is  clear  that 
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ionization,  recombination,  attachment,  and  multiply- 
charged  ions  may  be  ignored  for  the  systems  considered  here. 
This  does  not  imply  that  these  processes  are  unimportant, 
but  only  that  their  effects  can  be  included  in  a self- 
consistent  manner  for  a given  electron  density.  In  this 
case,  it  is  justifiable  to  regard  the  electron  number 
density  as  a parameter.  As  with  the  electric  field,  this 
parameter  is  conveniently  expressed  in  terms  of  the  fraction- 
al ionization,  ne/n0  . In  keeping  with  the  stable  discharge 
criterion,  the  range  of  ne/n0  is  restricted  to  10  ^ to 
10”5.  These  values  are  also  representative  of  those  found 
in  numerous  real  systems. 

Based  upon  the  discharge  parameters  specified  above, 
the  largest  fractional  excited-state  density 
consistent  with  stable  discharge  operation  is  approximately 
10"5,  and,  the  effects  of  superelastic  collisions  can  be 
neglected  (Ref  15).  To  strengthen  this  argument,  numerical 
calculations  have  been  performed  for  a variety  of  excited 
state  populations.  The  effect  on  the  equilibrium  electron 
distribution  is  found  to  be  quite  small  (Ref  20). 

The  next  assumption  considered  is  that  of  constant 
Maxwellian  distributions  for  the  neutrals  and  ions.  Actually, 
the  specific  distribution  of  these  species  is  unimportant 
since  it  does  not  enter  into  Eq  (22).  Furthermore,  since  it 
has  been  shown  that  the  fractional  ionization  is  a parameter 
and  the  excited-state  populations  can  be  ignored,  the  number 
densities  of  these  species  are  constant.  What  is  important 
is  the  energy  transfer  between  electrons  and  these  species. 
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Only  in  this  context  are  the  distribution  functions  for  the 
neutrals  and  ions  important.  The  same  can  be  said  for  the 
assumption  that  the  ions  and  neutrals  are  infinitely  massive, 
since  this  really  implies  a zero  kinetic  temperature.  There- 
fore, all  of  these  assumptions  are  justified  provided t 

(a)  the  energy  transfer  between  electrons  and  other 
species  is  negligible, 

(b)  the  kinetic  temperature  change  of  the  neutrals 
and  ions  is  negligible  during  the  time  the 
electron  distribution  reaches  steady-state. 

Both  of  these  conditions  are  reasonable  in  view  of  the  large 
mass  differential  between  electrons  and  other  species.  The 
neutral  particles  can  gain  energy  only  through  collisions 
with  the  electrons  and  ions.  Collisions  with  ions  have 
little  effect  because  the  neutrals  and  ions  have  comparable 
temperatures.  Collisions  with  electrons  are  inefficient 
because  of  the  large  mass  difference.  The  ions  can  transfer 
energy  with  electrons  or  neutrals  and  also  gain  energy  from 
the  field.  The  field  is  relatively  ineffective  in  heating 
the  ions  because  of  their  large  mass.  The  energy  transfer 
with  neutrals  and  electrons  has  little  effect  by  the  same 
argument  used  above. 

Even  though  there  is  a negligible  energy  transfer 
between  electrons  and  ions,  these  collisions  do  affect  the 
heating  rate  of  the  electrons.  This  effect  has  been 
neglected  by  assuming  l) eL<<  • Since  the  number  densi- 

ties of  the  various  species  have  been  determined,  these 
quantities  can  be  calculated.  From  Appendix  D, 
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Therefore 


■.-8 


*.  = 5.928  x 10-°  and  Vm  = n0  Rm. 


Vm  = 5.928  x 1011  sec-1 


(38) 


It  can  be  shown  that,  on  the  average  (Ref  23*258) 
, _ ‘K2nV/t  / / ktV'i  , 

3 ‘Ikt J (— ; 


where  ni  is  the  total  number  density  of  ions  (cm“-^)  and 
T is  the  kinetic  temperature  of  the  electrons.  Taking 
n^  = 10  cm”'5  and  KT  = 2ev  (liberal  estimates  for  the 

parameters  considered  here) 


= 6.237  1010  sec”1 


Therefore,  even  under  the  worst  conditions  considered  here, 
^ei  is  approximately  an  order  of  magnitude  less  than  V 

When  F(v,t)  was  expanded  in  spherical  harmonics,  it 
was  assumed  that  the  anistropic  term  was  small  compared  to 
the  isotropic  part  of  the  distribution  for  "small"  electric 
fields.  To  check  the  validity  of  this  assumption,  the 
electron  drift  velocity  was  calculated  numerically  for  all 
parameter  ranges  considered  here.  In  all  cases,  the  drift 
velocity  was  at  least  an  order  of  magnitude  less  than  the 
thermal  velocity  of  the  electrons. 

The  assumptions  of  a single  excitation  level,  constant 
momentum  transfer  rate,  and  constant  excitation  rate  are 
gross  approximations  to  reality.  These  assumptions  have 
been  made  to  simplify  the  Boltzmann  equation.  The  goal 
is  to  derive  an  approximate  analytic  solution  tc  the 
Boltzmann  equation  under  these  conditions.  Hopefully,  this 
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solution  will  predict  general  characteristics  of  the  system 
that  are  insensitive  to  the  detailed  cross  sections.  For 
numerical  solutions  to  the  Boltzmann  equation,  these  assump- 
tions are  not  necessary.  Another  problem  then  arises  in  that 
accurate  cross  sections  (especially  excitation  cross  sections) 
are  frequently  unavailable. 

Lastly,  the  assumptions  used  in  Appendix  A to  derive 
the  Fokker-Planck  equation  should  be  mentioned.  The  validity 
of  the  "weak-interaction"  approximation  is  established  in 
Appendix  A.  Correlation  effects  are  neglected  in  the  deriva- 
tion. Furthermore,  the  expansions  used  are  carried  out  to 
second  order  only  and  no  estimates  are  given  for  the  error 
introduced  by  higher  order  terms.  The  justification  of  these 
assumptions  is  non-trivial.  It  can  be  shown,  however,  that 
these  effects  are  quite  small.  A detailed  discussion  of 
correlation  effects  is  given  by  Kaufman  (Ref  7i295)»  Higher 
order  coefficients  are  estimated  by  Gasiorowicz,  et  al  (Ref  9). 

Summary 

The  Boltzmann  equation,  subject  to  a number  of 
assumptions,  has  been  cast  into  a relatively  simple  form. 

The  equation  includes  the  effects  of  a constant  electric 
field,  elastic  collisions,  electronic  excitation  and 
electron-electron  collisions.  The  equation  is  applicable 
for  e/a0  values  from  10_1®  to  10“*^  volts-cm^  and 
°e/f\ „ values  from  10“^  to  10“^.  The  relevant  equation 
is  Eq  (22),  with  the  Fokker-Planck  electron-electron  colli- 
sion term  given  by  Eqs  (33)  through  (37)  and  Eq  (31)* 


Ill  Approach 


Equation  (22)  is  a non-linear  integro-differential 
equation  for  the  steady-state  electron  distribution  function 
n(E).  The  non-local  term,  n(£+x),  further  complicates 
matters.  This  equation  cannot  be  solved  exactly  by  any 
known  techniques.  An  approximate  solution  is  possible  if 
the  electron-electron  collision  term  is  ignored.  To  extend 
the  solution  to  include  these  effects,  the  term  must  be 
approximated.  It  is  not  obvious  from  Eqs  (33)  through  (36) 
how  best  to  approximate  this  term.  Therefore,  the  general 
characteristics  of  the  Fokker-Planck  term  are  investigated 
first.  This  provides  direction  in  formulating  the  approxi- 
mate relaxation  term  and  also  an  indication  of  its  short- 
comings. Finally,  Eq  (22)  is  solved  approximately  using 
the  relaxation  terra  to  include  electron-electron  collisions. 

General  Characteristics  of  the  Fokker-Planck  Term 


From  physical  principles,  electron-electron  collisions 
are  expected  to  have  the  following  features  1 


(a)  The  total  electron  number  density  is  conserved. 

(b)  The  average  electron  energy  is  conserved. 

(c)  A Maxwellian  distribution  is  unaltered. 

It  is  straightforward  to  show  that  Eq  (33)  has  these 
features.  The  change  in  the  total  number  density,  due  to 
electron  collisions,  is  given  by 


it 
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where  Eq  (33)  has  been  used.  Since  n(£)  vanishes  at  zero 
and  infinity,  J(£)  also  vanishes  at  these  limits  by 
Eq  (3^).  Therefore,  the  Fokker-Planck  collision  term 
conserves  total  electron  number  density. 

Similarly,  the  change  in  the  average  energy  £ is 


_ jl  r 

at  “ J ' 
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Using  integration  by  parts  and  the  fact  that  J(£)  vanishes 
at  zero  and  infinity,  the  above  equation  becomes 


or  using  Eqs  (33)  through  (36) 
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The  second  integral  in  this  equation  can  be  integrated 
by  parts  as  follows i 

£ 

^ du-n  (€^£ 


€ 


Therefore 


The  first  integral  may  be  written  in  the  following  formt 


If  these  expressions  are  substituted  into  Eq  (39) » "the 

3 F 

right  side  vanishes  identically.  Therefore,  : O , and 
the  average  electron  energy  is  conserved  by  the  Fokker-Planck 


en 


The  last  physical  feature  to  be  verified  i6  station- 
arity  of  a Maxwellian  distribution  n(£),  where 


i/t  -£/kt 

* £ a 


(40) 


The  verification  is  straightforward  and  will  not  be  given  in 
detail.  Substituting  n(£)  into  Eqs  (34)  through  (36)  and 
integrating  the  first  integral  in  Eq  (35)  by  parts,  J(£) 
vanishes  for  all  £ . As  a result,  using  Eq  (33) 
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and  a Maxwellian  distribution  is  unaltered. 

Therefore,  the  Fokker-Planck  form  of  the  electron- 
electron  collision  term  retains  the  three  main  physical 
features  of  these  interactions.  Parenthetically,  this  gives 
further  credibility  to  the  assumptions  used  to  derive  the 
Fokker-Planck  equation. 

An  inspection  of  Eqs  (35)  and  (36)  reveals  that  the 
coefficients,  P and  Q,  involve  integrals  of  n(£).  This 
not  only  makes  the  equations  non-linear,  but  also  introduces 
a non-local  effect.  Therefore,  |ee  depends  not  only 

on  n or  its  derivatives  at  a given  energy,  but  also  on  the 
shape  of  the  distribution  function  over  the  entire  energy 
range.  This  is  reasonable  (and  necessary)  since  electron- 
electron  collisions  cannot  change  the  number  density  or 
average  energy.  Several  techniques  have  been  used  to  approxi- 
mate this  effect  (Ref  17). 
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Unfortunately,  attempts  to  include  the  non-local 
effects  in  a simple  collision  term  failed.  To  retain  the 
physical  features  discussed  above,  a Krook-type  relaxation 
term  is  assumed i 

^n(€)|  _ _ w(Q- nCO 

d-t  let  tret') 

where  is  the  characteristic  time  in  which  the  distri- 

bution relaxes  to  a Maxwellian  as  a result  of  electron- 
electron  collisions. 

In  general,  X.  must  be  considered  a function  of 
energy.  There  is  no  reason  to  expect  the  distribution  to 
relax  uniformly  at  the  same  rate.  In  fact,  numerical  solu- 
tions show  that  the  high  energy  portion  of  the  distribution 
relaxes  much  slower  than  the  low  energy  portion  (Ref  16). 

A qualitative  description  is  provided  by  considering  similar 
distributions  with  different  average  energies.  The  result 
ie  f (Ref  6*120) . It  is  not  clear  that  this 

relationship  is  generally  valid  for  a given  distribution. 
However,  this  form  for  ■£'<£')  does  insure  that  the  relaxa- 
tion rate  decreases  with  increasing  energy. 

Attempts  to  include  an  energy-dependent  relaxation 
rate  and  retain  all  three  physical  features  have  also  been 
unsuccessful.  To  insure  the  relaxation  term  conserves  number 
density  and  average  energy,  n(fc)  must  be  altered  in  a 
manner  determined  by  t' (£) . With  this  alteration,  n(£) 
no  longer  retains  its  true  Maxwellian  nature  (with  the  same 
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number  density  and  average  energy  as  the  actual  distribution). 
The  basic  problem  is  that  for  an  energy-dependent  relaxation 
rate,  the  physical  features  overspecify  n(£).  Therefore, 
to  employ  the  simple  relaxation  term  a constant  't'  is 
used , 

The  value  of  /t/  used  here  is  the  "self-collision 
time"  t'ce  derived  by  Spitzer  (Ref  24 1 1 33 ) ; it  is  basically 
the  average  time  it  takes  an  electron  with  average  energy  to 
undergo  a 100ft  change  in  energy  and  suffer  a 90°  deflection 
(Ref  16).  The  self-collision  time  may  be  expressed  as 


o.zfcfc  T 3/tm 
ne  In  A 


(sec) 


(41) 


where  T is  the  kinetic  temperature  of  the  electrons  (°K). 
The  final  form  for  the  relaxation  term  is 


3rv(Ol  _ _ n_a)  - (42) 

dt  lee 

It  is  easily  verified  that  this  form  of  the  electron- 
electron  collision  term  retains  all  the  physical  features 
discussed  earlier,  provided  n(£)  has  the  same  number 
density  and  average  energy  as  n(g).  It  is  anticipated  that 
this  term  will  produce  a higher  than  normal  relaxation  rate 
at  high  energies  and  a lower  than  normal  rate  at  low  energies. 

A problem  does  arise  when  using  Eq  (42)  for  finite- 
differenced  numerical  solutions.  In  this  case,  n(£)  is 
not  well  defined  by  Eq  (40).  The  method  used  to  circumvent 
this  problem  is  given  in  Appendix  C. 
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Approximate  Solution  of  the  Boltzmann  Equation 


The  Boltzmann  equation  may  be  solved  approximately 
if  the  relaxation  term  is  used  for  electron-electron  colli- 
sions. Combining  Eqs  (22)  and  (42),  the  Boltzmann  equation 
becomes 

■j  i - £ Xl]  - n0R(£.+  x^A(£4-*^ 

-ruR(0n(O-  (43) 

A fundamental  problem  arises  immediately  in  that  n(£) 
depends  on  the  steady-state  average  energy,  which  is  unknown 
at  this  point.  This  problem  is  not  peculiar  to  the  relaxation 
term;  the  Fokker-Planck  term  also  depends  on  the  average 
energy  through  the  quantity  UA  . However,  electron-electron 
collisions  usually  have  little  effect  on  quantities  averaged 
over  the  entire  distribution  function  (Ref  1»22).  Therefore, 
Eq  (43)  could  be  solved  using  the  techniques  mentioned 
earlier  and  ignoring  electron-electron  collisions.  The 
average  energy  obtained  from  this  solution  could  then  be 
used  to  evaluate  n(£). 

The  solution  to  Eq  (43)  will  be  restricted  to  the 
region  X . In  this  region  it  is  reasonable  to  assume 
r'oR(t+x')n(€+x)<<n0R(£)r)(€).  The  numerical  solutions  presented 
in  the  next  section  show  that  this  assumption  is  good  until 
electron-electron  collisions  completely  dominate  inelastic 
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collisions.  The  equation  may  be  cast  into  a more  convenient 
form  by  the  following  change  of  variable i 


In  addition 


= -It' 


h(s')  r 


ne  £‘l2  n (O 


a(0  = 


c 1,1  ' 

1 e 


where 


Jh.  l3/2. 

o-  -m  b 
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(OciE  = 


With  these  changes  and  assumptions,  Eq  (43)  takes  the 
following  form« 


- A,  K = ^ 


where 


A = [\  (^oRo  + ^tj] 


(45) 


- 


24  a. 


-4-b/V 


(46) 


t 1, 


ce 


Equation  (44)  can  be  solved  using  a combination  of 
the  WKB  method  and  variation  of  parameters.  Unfortunately, 
the  complete  solution  involves  integrals  which  cannot  be 
evaluated  in  closed  form.  Since  this  approach  could  be 
useful  for  future  work;  the  results,  though  incomplete,  are 
included  in  Appendix  E. 

An  approximate  solution  can  be  obtained  by  ignoring 
the  second  derivative  in  Eq  (44).  In  this  approximation, 
n is  given  by 


. _ 2.  4 o- 

TT^r  e 
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ce 


Using  Eq  (45),  the  solution  may  be  expressed  as 


CL 


l + n.fi.t 


ce 
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Equation  (47)  represents  a self-consistent  solution  provided 


di 

_n 

U 


A 1 


« 
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(48) 
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This  requirement  follows  directly  from  Eq  (44) . Differen- 
tiating Eq  (4?)  and  substituting  the  result  into  Eq  (48) 
leads  to  the  following  requirement: 


In  the  next  section,  three  numerical  solutions  to 
Eq.  (43)  are  compared  with  solutions  using  the  Fokker-Planck 
collision  term.  It  is  easily  verified  that  the  above 
requirement  is  satisfied  for  these  cases.  Transforming 
Eq  (47)  back  into  energy  space  yields  the  following  approxi- 
mate solution  to  Eq  (43)  for  ElXt 

'/*■  -£/kt 


h(0  = 


ne  £ 


I +n0R0Tte 


(49) 


Since  Eq  (49)  involves  no  arbitrary  constants,  the  complete 
solution  to  Eq  (43)  (subject  to  the  imposed  physical 
restraints)  will  generally  be  discontinuous  at  the  excita- 
tion energy.  This  behavior  and  several  other  features  are 
investigated  in  the  next  section.  Since  a discontinuous 
distribution  function  is  not  physically  acceptable,  no 
attempt  is  made  to  extend  the  solution  below  the  excitation 
energy. 

Summary 

A simple  relaxation  term  has  been  formulated  that 
retains  the  three  prominent  physical  feature  of  electron- 
electron  collisions.  Non-local  effects  have  been  disre- 
garded. Using  this  term  for  electron-electron  collisions, 


I 


33 


an  approximate  solution  to  the  Boltzmann  equation  is 
obtained  for  energies  greater  than  the  excitation  energy. 
The  relaxation  term  is  given  by  Eq  (3*0  » the  approximate 
solution  is  given  by  Eq  (49). 


34 


IV  Analysis  and  Comparison 


In  this  section,  the  relaxation  form  of  the  electron- 
electron  collision  term  is  compared  to  the  Fokker-Planck 
collision  term.  The  temporal  characteristics  of  the  two 
terms  are  analyzed  first.  Next,  steady-state  solutions  of 
the  Boltzmann  equation  are  analyzed  for  three  sets  of  £/n0 
and  rWn©  values.  The  Fokker-Planck  and  relaxation 
results  are  computed  numerically  for  all  cases.  These 
computations  employ  a finite-differenced  energy  axis  with 
a 20  ev  maximum  energy  and  .5  ev  energy  interval.  The 
analytic  results  are  calculated  using  Eq  (49)  for  the 
three  steady-state  comparisons. 

Temporal  Characteristics 

The  initial  distribution  is  chosen  to  be  a Gaussian 
with  a 10  ev  average  energy  and  1 ev  standard  deviation 
(Fig.  1).  The  temporal  characteristics  of  the  two  terms 
are  compared  in  Figs.  2 through  5»  Only  the  effects  of 
electron-electron  collisions  are  included  in  these  calcu- 
lations. The  numerical  procedure  described  in  Appendix  B 
is  used  to  calculate  the  Fokker-Planck  term?  the  relaxation 
term  is  calculated  using  the  procedure  in  Appendix  C. 

As  expected,  the  relaxation  term  gives  a higher  than 
normal  rate  at  high  energies  (Fig.  2).  However,  the  rate  is 
also  higher  at  low  energies.  This  is  just  the  opposite  of 
the  predictions  in  Section  III,  This  result  can  be  explained 
by  considering  the  Maxwellian  used  in  the  relaxation  calcu- 
lation, It  has  been  verified  numerically  that  this 
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Maxwellian  (time-independent  in  this  case)  corresponds 
exactly  to  the  final  steady-state  distribution  of  the 
Fokker-Planck  term.  This  distribution,  shown  in  Fig.  5> 
is  essentially  a horizontal  line  for  the  log  scale  used 
here.  Since  the  initial  Gaussian  is  symmetrical  about  the 
average  energy,  the  relaxation  term  predicts  approximately 
equal  rates  for  low  and  high  energies. 

Therefore,  no  generally  valid  relationship  exists 
between  the  Fokker-Planck  rates  and  the  relaxation  term 
rates.  Furthermore,  an  energy-dependent  could  elimi- 

V v 

nate  this  problem  for  one  specific  distribution  only.  The 
non-local  effects,  neglected  in  formulating  the  relaxation 
term,  must  be  included  to  completely  eliminate  the  problem. 

The  form  of  the  distribution  at  t = 10~10sec  is 
shown  in  Fig.  3.  The  relaxation  term  exhibits  the  same 
features  as  before}  however,  the  diffusion  characteristics 
of  the  Fokker-Planck  term  are  now  evident.  At  this  point, 
the  distribution  has  spread  (slightly  weighted  toward  the 
higher  energies)  about  the  average  energy  and  the  number  of 
electrons  with  the  average  energy  has  decreased  slightly. 

Figure  4 shows  the  form  of  the  distribution  at 
t = 10”^sec.  This  time  closely  approximates  "tf  for  the 
parameters  assumed  here.  Obviously,  the  Fokker-Planck  and 
relaxation  distributions  are  considerably  different.  However, 
Fig.  4 does  indicat'  that  in  either  case,  'C  is  a reasona- 
ble  time  scale  for  electron-electron  collisions  to  relax  a 
distribution  to  a Maxwellian. 
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The  steady-state  condition  is  shown  in  Fig.  5*  No 
later  times  are  shown  because  the  numerical  integration 
procedure  was  unable  to  advance  the  relaxation  term  any 
further.  It  appears  that  the  relaxation  term  is  not  well 
suited  to  numerical  integration.  As  mentioned  previously, 
the  Maxwellian  used  in  the  relaxation  calculation  is  identi- 
cal to  the  Fokker-Planck  steady-state  distribution.  There- 
fore, both  approaches  should  yield  identical  steady-state 
results.  This  problem  was  not  encountered  when  the  initial 
distribution  differed  only  slightly  from  a Maxwellian.  In 
this  case,  both  methods  result  in  identical  steady-state 
distributions.  The  large  average  energy,  coarse  energy 
zoning  and  relatively  small  maximum  energy  could  be  a 
factor.  Time  did  not  permit  a closer  investigation  of  this 
problem. 

Steady-State  Comparisons 

Steady-state  distributions  for  three  sets  of  E/ne 
and  n*/n0  values  are  shown  in  Figs.  6 through  8.  The 
Fokker-Planck  and  relaxation  curves  represent  numerical 
calculations.  The  analytic  curve  represents  the  evaluation 
of  Eq  (49)  for  the  appropriate  parameters.  The  excitation 
threshold  energy  is  8 ev  in  all  cases.  For  completeness, 
the  results  of  the  numerical  calculations  are  tabulated  in 
Appendix  F. 

There  are  two  characteristic  features  of  all  three 
cases.  First,  the  Fokker-Planck  and  relaxation  approaches 
agree  well  at  energies  below  the  excitation  threshold.  This 
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is  to  be  expected  because  there  are  no  energy  loss  mechanisms 
in  this  range;  energy  gained  from  the  field  must  be  redistri- 
buted by  electron-electron  collisions.  Inelastic  and  super- 
elastic processes  are  treated  identically  in  both  numerical 
procedures.  Therefore,  the  precise  rate  of  electron-electron 
collisions  should  have  little  effect  on  the  steady-state 
solutions  for  energies  below  the  excitation  threshold. 
Secondly,  the  distributions  are  essentially  Maxwellian  over 
the  entire  energy  range.  The  Fokker-Planck  distributions 
exhibit  a lower  temperature  above  threshold.  The  relaxation 
distributions  can  be  characterized  by  a single  temperature 
over  the  v/hole  energy  range.  This  important  and  discouraging 
feature  will  now  be  investigated  more  thoroughly. 

The  effect  of  the  large  relaxation  rate,  character- 
istic of  the  relaxation  term  at  high  energies,  is  evident 
in  all  three  cases.  In  fact,  the  inelastic  collision 
process  is  completely  overwhelmed.  The  only  departure  from 
a true  Maxwellian  is  a discontinuity  in  the  neighborhood  of 
the  excitation  threshold.  The  largest  discontinuity  occurs 
in  Fig.  6,  This  case  corresponds  to  the  smallest  field  and 
fractional  ionization  considered.  The  smallest  discontinuity 
occurs  in  Fig.  7,  This  case  corresponds  to  a small  field 
and  large  fractional  ionization.  Therefore,  the  approaches 
agree  better  as  the  number  density  increases  and  electron- 
electron  collisions  become  more  important. 

These  results  may  be  explained  as  follows.  First,  the 
relaxation  term  tends  to  overemphasize  electron-electron 
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collisions  in  all  cases.  Since  there  is  an  upper  limit  to 
this  effect  (a  true  Maxwellian  distribution),  the  relaxation 
theory  and  Fokker-Planck  theory  agree  better  at  higher 
fractional  ionizations.  Secondly,  the  relaxation  theory 
has  been  forced  to  conserve  number  density  and  average 
energy.  This  requirement,  together  with  the  rapid  relaxa- 
tion rate  at  low  and  high  energies,  produces  a discontinuity 
in  the  neighborhood  of  the  excitation  threshold.  Therefore, 
the  magnitude  of  this  discontinuity  can  decrease  as  the 
actual  distribution  approaches  a Maxwellian. 

The  average  energies  of  the  distributions  are  given 
in  Tables  III  through  VIII.  These  values  indicated  that 
the  two  approaches  agree  to  within  five  per  cent.  The 
largest  difference  occurs  for  ne/n0  = 10”^  and  Vn0  = 10 
These  are  also  the  conditions  for  which  the  low  energy  portion 
of  the  distribution  is  affected  most.  Nevertheless,  the 
assumption  made  concerning  averages  over  the  entire  distri- 
bution function  seems  to  be  justified. 

The  analytic  solution  derived  in  Section  III  is  also 
plotted  in  Fig.  6 through  8.  Obviously,  the  agreement 
between  the  analytic  and  numerical  solutions  is  quite  good. 

In  addition,  it  can  be  concluded  that  the  problems  encountered 
with  numerical  solutions  using  the  relaxation  model  are  not 
numerical  in  origin  but  result  from  the  discontinuity  at 
the  excitation  threshold.  It  is  surprising  that  the  numeri- 
cal integration  procedure  is  able  to  handle  the  problem  as 
well  as  it  does.  The  slight  difference  in  slope  between 
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the  numerical  and  analytic  solutions  could  result  from 
neglecting  the  second  derivative  in  Eq  (44).  Another 
possible  explanation  is  that  the  approximation  made  in 
solving  this  equation  is  not  good  for  very  large  energies. 

Furthermore,  Eq  (49)  reveals  that  changing  the 
magnitude  of  V would  not  alter  the  form  of  the  relaxa- 

L c 

tion  solution;  the  magnitude  of  the  discontinuity  would 
be  affected  but  the  overall  form  would  be  unchanged.  This 
argument  is  valid  provided  the  assumptions  made  in  deriving 
Eq  (49)  are  justified.  Since  an  energy-dependent 
would  change  the  nature  of  the  analytic  solution,  it  is 
possible  that  such  an  approach  would  lead  to  better  agree- 
ment with  the  Fokker-Planck  theory. 
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V Conclusion 


The  analysis  and  comparisons  presented  in  Section 
IV  lead  to  several  logical  conclusions  concerning  the 
relaxation  approach  to  electron-electron  collisions.  The 
approach,  as  presented  here,  is  not  applicable  for  the 
electric  fields  and  fractional  ionizations  commonly  used 
in  discharge  lasers.  The  high  energy  portion  of  the 
distribution  function  is  of  primary  importance  in  these 
devices.  The  relaxation  collision  term  yields  distribu- 
tions too  large  by  several  orders  of  magnitude  in  this 
range , 

The  three  most  serious  defects  in  the  relaxation 
approach  are  the  following i 

(a)  The  high  relaxation  rate  at  high  energies 
completely  dominates  the  more  important 
inelastic  processes. 

(b)  Energy  diffusion,  an  important  characteristic, 
is  ignored. 

(c)  The  steady-state  distribution  is  in  general 
discontinuous  at  the  excitation  threshold. 

The  analytic  solution  using  the  relaxation  approach 
indicates  that  a different  relaxation  rate  would  not 
significantly  improve  the  situation.  However,  an  energy- 
dependent  rate  could  possibly  lead  to  better  agreement 
with  the  Fokker-Planck  theory. 

Regardless  of  the  approach,  diffusion  effects 
appear  to  play  a key  role  in  representing  the  physical 
situation.  It  is  not  clear  how  to  include  this  feature 


in  a collision  term  with  reasonable  simplicity.  Finally, 
a simple  but  accurate,  numerical  procedure  might  be  more 
feasible  and  valuable  than  an  approach  based  on  drastic 
assumptions. 
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Appendix  A 


Derivation  of  the  Fokker-Planck  Equation 

The  purpose  of  this  appendix  is  to  derive  an  expression 
that  adequately  describes  the  effects  of  Coulomb  collisions 
on  the  electron  velocity  distribution  function.  The  Boltzmann 
collision  term  is  based  on  the  assumption  that  particle 
interactions  are  short-range,  binary  encounters.  For  electron- 
neutral  interactions,  this  assumption  is  usually  justified. 

To  describe  Coulomb  collisions,  the  effect  of  simultaneous 
interactions  between  many  particles  must  be  considered. 

It  is  not  possible  to  solve  this  many-body  problem 
exactly,  taking  all  the  particle  correlation  effects  into 
account.  However,  if  the  interactions  are  assumed  "weak", 
these  simultaneous  interactions  may  be  regarded  as  the  sum 
of  many  uncorrelated  binary  collisions.  The  following 
considerations  clarify  and  lend  credibility  to  this  approach. 

A weak  interaction  is  one  in  which  the  momentum  change 
of  a particle  is  small  compared  to  the  total  particle  momen- 
tum. In  highly  ionized  plasmas,  where  Coulomb  collisions 
must  be  considered,  the  cumulative  effect  of  many  weak 
interactions  is  more  effective  in  deflecting  an  electron 
than  a single  strong  interaction.  It  can  be  shown  (Ref  l4« 

294)  that  weak  interactions  dominate  strong  interactions  by 
a factor  of  8lnA  , where  JY  is  the  Coulomb  factor  defined 
in  Section  II.  Since  In  A usually  lies  between  10  and  20, 
weak  interactions  represent  the  overall  Coulomb  collisional 
effect  to  several  orders  of  magnitude. 
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The  simultaneous  interactions  may  now  be  regarded  as 
simultaneous  small  perturbations  and  the  total  effect  obtained 
by  adding  the  individual  perturbations  linearly.  Using  these 
ideas,  the  form  of  the  Coulomb  collision  term  can  be  derived 
by  several  independent  methods.  The  most  general  procedure 
starts  with  the  Liouville  equation  (Ref  11:115).  Another 
method  involves  expanding  the  Boltzmann  collision  integral  in 
powers  of  the  momentum  interchange  (Ref  8).  The  method  used 
here,  similar  to  that  used  in  the  study  of  Brownian  motion 
(Ref  23:31),  results  in  the  Fokker-Planck  equation.  This 
equation  correctly  describes,  to  second  order,  the  temporal 
evolution  of  the  electron  velocity  distribution  function 
resulting  from  many  weak,  uncorrelated  binary  encounters. 

Let  F(v,t)  represent  the  spatially  homogeneous  electron 
velocity  distribution  function.  TCy  - Ay,  A y)d  (Ay)  is  defined 
as  the  probability  that  an  electron  changes  its  velocity 
from  \/-AY  to  V in  a time  At  as  a result  of  Coulomb 
interactions.  Since  H^( V-Ay,A^)cl(Ay)  is  a probability, 
the  following  relation  must  hold* 


^ H^V-Ay,  Av)d(Ay)  - | (50) 

The  distribution  function  at  time  t +At  is  related  to  the 
distribution  at  time  t by 


F(v,  -t  + Al)=  \ FCV'A^l)  T(y-Ay)Av)4(AV)  (51) 


k 
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Since  weak  interactions  are  being  considered,  A v is  small 
in  a small  time  interval  At.  The  terms  in  Eq  (51)  can  be 
expanded  as  follows* 


F(v>t  + At)  = FU,-t)+  At 


F(v-Av)t')  = F(v)t')  - Z 

+ ^AV;fiVi®r+ 

ntix.flit)  =V(v,w) 


(52) 


(53) 
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where  Av/^  is  the  cth  component  of  A v.  If  Eqs  (52) 
through  (54)  are  substituted  into  Eq  (51),  the  result  is 

F(v.O  + 

-|y.  ^AVt-  t’C'ijAv)  d(Av') 

-Z  1^-  | A \ZC  Y(v,  AK)d(dv) 

+ T XI  F(v,0  Pa  V{.  A V;  T(v',AV)d(Av^ 

4.  _L  V"  p 

2.  Z-  Avc  AVj  ^(Y,  AV)d(Ay')  (55) 
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The  first  terms  on  the  left  and  right  sides  of  Eq  (55) 
cancel  identically  by  Eq  (50)*  Defining 


~ ci(A 


<AV-A V,->  = AV'jfCy.Av^CAv) 


Eq.  (55)  may  be  rewritten  as 


^X^=-I!7rFCv.O<*v;> 

dt  ; 6 


R^<^.iVi)  (58) 

ii5  t 


Equation  (58)  is  the  Fokker-Flanck  equation,  and  is 
phenomenological  until  explicit  expressions  are  derived 
for  Eqs  (56)  and  (57).  The  form  of  these  equations  is 
considered  in  Section  II  and  will  not  be  derived  here. 

Short  of  this,  several  general  remarks  can  be  made  concerning 
the  Fokker-Planck  equation. 

Equation  (58)  is  basically  a conservation  equation  in 
velocity  space.  The  terms  <AVC>  and  <AVtAV->  may  be 
interpreted  as  diffusion  coefficients  (Ref  24:125).  It  is 
important  to  note  that  these  coefficients  are  time-ensemble 
average  quantities  and  differ  dimensionally  from  the  corres- 
ponding ensemble  average  quantities.  The  quantity 
has  been  termed  the  coefficient  of  dynamical  friction 
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(Ref  2).  It  represents  the  average  deceleration  experienced 
by  an  electron  of  velocity  v due  to  Coulomb  interactions 
(Ref  8).  Similarly,  the  quantity  <AV,;avj-^  is  called  the 
coefficient  of  diffusion  tensor.  This  term  describes  the 
random  fluctuations  about  the  average  velocity  (Ref  23*249). 

In  conclusion,  it  should  be  mentioned  that  applica- 
tions of  the  Fokker-Planck  equation  are  not  restricted  to 
Coulomb  interactions  in  plasmas.  In  fact,  the  theory  has 
been  extensively  used  to  analyze  problems  in  stellar 
dynamics  (Ref  2,3,4).  A mathematically  rigorous  discussion 
is  given  by  Haken  (Ref  10). 


Appendix  B 

Numerical  Solution  of  the  Fokker-Planck 
Electron-Electron  Collision  Term 
(Subroutine  Funse ) 

A significant  part  of  this  thesis  effort  has  been  the 
development  of  a computer  code  to  solve  the  Boltzmann  equa- 
tion under  a variety  of  conditions.  This  work,  in  coopera- 
tion with  three  other  students,  resulted  in  a computer 
program  called  NG3 . The  operating  procedure  for  NGB  is  to 
be  documented  in  an  AFIT  Technical  Report.  By  mutual  agree- 
ment, the  details  of  the  various  segments  of  NGB  are  presented 
as  Appendices  in  each  individual's  thesis.  This  appendix 
fulfills  this  requirement  for  the  electron-electron  collision 
term. 

The  equations  used  to  numerically  evaluate  the  electron- 
electron  collision  term  are  derived  and  presented  in  a form 
suitable  for  direct  computer  coding.  Particular  attention  is 
given  to  units  and  numerical  factors.  The  method  used  here 
was  developed  by  Proctor  and  Canavan  (Ref  19) , although  the 
notation  parallels  that  of  Rockwood  (Ref  20).  Subroutine 
Funee,  the  coding  used  fox’  electron-electron  collisions  in 
NGB,  is  inclxided  at  the  end  of  this  appendix. 

From  Section  II,  the  electron-electron  collision  term 
may  be  expressed  as 
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where  n(£,t)d£  is  the  number  density  of  electrons  with 
energies  between  £ and  £+d£  at  time  t,  and 


J=*[p(zV 

i OC 

p e Z £ ,ll~  J^VC  r»  C x, -L')  d x + Z £ J x~  /z  rt  (*,*)  d 

6 t 

rc 

Q - 3 1 J n£x,t)  4 x 


(34) 

(35) 


(36) 


Z.  u / 2.  \'/z  j 

^ - 3 7T  e (m  ) | n A 


A - -3,  / KaT?y/a. 
2e3  ktt  hq  J 


(37) 

(31) 


and  ne  = total  number  density  of  electrons 


Units  and  numerical  factors  will  be  considered  first. 
Energies  are  assumed  to  be  in  electron  volts  (ev*s)  and 
number  densities  per  cubic  centimeter  (cirT^) , The  other 
quantities  of  interest  must  have  the  indicated  units  to 
insure  consistency! 


n(Z,±)  ^ cm'3-  ev' 


2 t 


- 3 - I 


^ c m"  e.v  _ s - * " 1 


sec 


ee. 


- 3 3/?  I 

d,  cm  -c.v  -sec1 


A - dimension  less 
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From  Eq  (31 ) , A can  be  rendered  dimensionless  by  entering 
the  electron  charge  in  esu’s,  ng  in  cm  J and  KT  in 


ergs.  Denoting  the  average  energy  by  £ , the  result  is 

. . 97  r Vz  (er9s)v 

A = 4.15?64  x 102' 


Since  1 ev  = 1.60210  x 10~*2  ergs 

A ■ 


and 


InA  = 22.85519  + In  { -^7!  J 


(59) 


From  Eq  (37) » with  the  electron  mass  entered  in  grams  and  the 

3 3/2  -1 

charge  m esu’s,  the  units  for  are  cm'/-ergs'v  -sec 
After  converting  from  ergs  to  ev's,  the  following  result  is 
obtained  1 


* = (2.57540  x 10~6)(nA  (cm^-ev^-sec"1) 


(60) 


Equation  (33)  represents  only  one  term  in  the  Boltzmann 
equation.  Therefore,  the  numerical  procedure  used  to  evaluate 
this  term  must  be  consistent  with  the  method  used  to  solve 
the  Boltzmann  equation.  A brief  description  of  the  method 
used  in  NGB  is  in  order. 

The  Boltzmann  equation  is  finite-differenced  into  a 
set  of  N discrete  energy  intervals  of  width  A£  , The 
resulting  time  dependent  equation  is  then  integrated  using 
a forward-marching  sixth  order  Gear- type  algorithm.  The 
initial  distribution  is  Maxwellian  with  reasonable,  but 
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arbitrary,  parameters.  For  further  details,  the  Technical 
Report  describing  NGB  should  be  consulted. 

Equation  (33)  must  now  be  expressed  in  a form  consis- 
tent with  the  above  differencing  scheme,  A diagram  of  the 
energy  axis  is  given  in  Fig.  9.  J.  represents  the  electron 
number  density  flux  (along  the  energy  axis)  resulting  from 
electron-electron  collisions. 

The  differenced  forms  of  Eqs  (33)  through  (36)  are 
given  below: 


t ~ 

“ A £ 

(61) 

Ji* 

' ntui  + ift 

"ft*, -"ft') 

-^lA  a / 

^ vr 

A £ / 

- 

(62) 

1 1 1 

l A £ + 2.  £ ^ 

(63) 

Qft 

i - I 

A £ 

(64) 
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To  insure  a constant  total  electron  number  density,  the 
following  relations  are  imposed* 


j:  =r*  = o 

j;  = jl 


Substituting  Eq  (62)  into  Eq  (66)  yields 


hfct-  nfe-!  _ r\  g - n fe-i 
‘tefe-vx  ae  . 


Q^+-Q  &_(\  / n fe  + *"*  ^ ~ 1 


Equations  (62)  and  (6?)  are  now  substituted  into  Eq  (6l) 


to  give 
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where  Eq  (65)  has  been  used  for  k=l  and  k=N . If  this 
set  of  equations  is  expanded  and  the  coefficients  of  the 
various  n^  collected,  the  following  set  of  equations  is 
obtained  1 


ft,  ~ _2dt 


{[CP^RXirh;^)-  -r~ 


- tit *•  ae~) - a**aa~3  oft. 

-(pfeo-  n* 

4(Pv^uX^p:/t-  h)  - "kJ 


•N  2A  £ 


1 {[(?„- UC?tg,. 


J_\  Qei  + Qv-i 

/>£  J ?. 


This  may  be  written  more  compactly  by  defining  the  following 
quantities  s 


a6  2 AC  {(  ^V5"  Ptc*.X  4 £fetv2  4 A t)  " 


i_\  % + Oi v 


bi  = 

K 


“2Ae (( + Fk-.Xnefc.*/* as)" 


Using  Eqs  (69),  Eqs  (68)  can  be  expressed  in  the  folio-wing 
simple  form: 


* 1 ' / 

rt,  - - Q_,  n ( 


H - aVi  " (ak  4*  nl  v nfL 


^ N- 1 " t#N  *N 


(70) 


Equation  (70)  is  used  to  evaluate  nk  in  Funee.  It  remains 
to  find  explicit  expressions  for  a ^ and  bj  . The  following 
definitions  considerably  simplify  the  derivation: 


4=0 


I + o.  25- 


A £ 


(71) 
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(72) 


Using  Eq  (71),  Eqs  (69)  may  be  written  in  the  following  formi 


Of 
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f •+  Ok  4‘  Q&+1I 

(•VPl £ j 
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<x 


k ~ 2^£ 


lCPt  - PJ 


+ 


] 


(73) 


Using  Eq  (72).  Eqs  (6 3)  and  (64)  become 


Pt  = ti'XO-H,.  , /> 

JL-  i i = i K-l}V~ 


N 

0^=3Efe.',/iA£  £L  nLH| 

£=  i 


(74) 


Equations  (73)  can  now  be  written  in  the  following  convenient 
form  t 


H--  !>{&%  h.**  *x  Hblil 


+tiHii.'hti-H£,i)]‘ut-0-7S(ekk  HU+  H 6*1, n)]  nt  ( 75) 


Defining  two  matrices  A „ and  B_  . , Eqs  (75)  and  (?6) 

ft  fl»x. 

become 


~ Z Ag  8 n_ 


X-i 


lfc(4_r,JL 


^ = Z_Br,  nhjL 

i=l 


(77) 
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+[£*0-Hfcl>t)+ElM(i-H*Je;%;}  (78> 


The  corresponding  relation  for  3„  „ is  not  needed 
explicitly.  This  will  become  clear  in  what  follows.  If 
Eqs  (77)  are  inserted  into  Eqs  (70)  and  summed  over  all  k, 
the  resulting  expression  vanishes  for  any  ^ and  ^ 
whatsoever.  Therefore 


and  the  total  electron  number  density  remains  constant. 
This  important  feature  is  retained  by  the  differenced 
expressions,  and  is  basically  a result  of  the  conditions 
imposed  by  Eqs  (65)  and  (66) . 


Unfortunately , when  Eqs  (70)  are  multiplied  by 
and  summed  over  all  k (with  Eqs  (77)  inserted  as  before) 
the  result  is 


the  average  energy  remains  constant  only  if  the  right  side 
of  Eq  (79)  vanishes.  To  insure  that  this  feature  is 
retained  by  the  differenced  expressions,  the  following 
conditions  must  be  imposed: 


B 


j.*- 


(80) 


B,,x=  s l <81> 

f\  , — O -for  all  S.  (82) 

IV  | a- 

Combining  Eq  (80)  and  Eq  (81)  leads  to  the  additional 
requirement  t 

A ZlO-Porall  Si  (8.3) 

nA,l 


The  differenced  expressions  may  now  be  reformulated  in  terms 
of  A»  ,,  a'r  and  b& . For  convenience,  the  results  so  far 
are  summeriz,ed  below: 
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and 


+foo-Htl>€^o-H,ixVl|  ™ 


It  was  shown  in  Section  III  that 


^ I - O 

^•t  tee. 


when  n(£,t)  is  Maxwellian.  The  differenced  expressions 
developed  thus  far  do  not  retain  this  property.  Consider 


a Maxwellian  distribution  n*  where 


n*  * ^ e"  fE/l<T 


and  K is  Boltzmann’ s constant.  It  is  easily  verified  that 

for  all  i and  -k.  From  Eqs  (?0)  and  (84) 
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The  first  sum  in  Eq  (86)  may  be  rewritten  using  Eq  (85) i 


A- I A- I 


tJ-i 


- y a t ^ ^ Vi-  - 

nn<-,n' 


X-  1 


i,e^j  nt+« 


;fA  C t,  £,.,V'*_  _ 
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Since  A11  = 0,  the  second  sum  in  Eq  (86)  can  be  rewritten 
also : 

N N 

ZA.iSA-  y^\, 

A- « XTT, 

It  follows  that  n,  vanishes  if: 


The  same  procedure  nay  be  applied  to  Eqs  (87)  and  (88). 
general  expression,  obtained  from  Eq  (8?),  is 


(89) 

The 


Equation  (90)  is  sufficient  to  guarantee  that,  when  n^  is 
J-iaxwelli an,  n^  = 0 for  all  k.  It  is  evident  from  Eq  (90) 

that  the  A-  . are  not  independent.  In  fact,  only  those 

K»x. 

above  the  main  diagonal  need  be  calculated  from  Eq  (78); 
the  remaining  terms  are  given  by  Eq  (90).  It  should  be 
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mentioned  that,  while  the  preceeding  arguments  have  been 
strictly  mathematical,  identical  results  can  be  obtained 
from  more  physical  arguments  (Ref  19). 

In  summary,  given  a set  of  n^,  n^  due  to  electron- 
electron  collisions  is  calculated  using  the  following 
equations : 

n,  = b*  rij_  - a.',  n,  - 1 

= ‘H-.^-rCoV  bUlnu. 

nN~  aN-i  hN-i  ” bw  (?o) 


where 


(84) 


(78) 


for  X ">  ft 
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V / 


(90) 
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This  procedure  strictly  conserves  total  electron 
number  density  and  average  energy.  For  a Maxwellian, 
n^  = 0 for  all  k.  Since  the  A matrix  depends  only  on 
the  differencing  scheme  and  not  on  the  distribution  function, 
A is  calculated  only  once  and  stored  for  subsequent  calcula- 
tions. As  a result,  execution  time  is  minimal. 

There  are  several  disadvantages  to  this  procedure. 

The  size  of  the  A matrix  is  given  by  the  square  of  the 
number  of  energy  intervals  used.  Thus,  storage  requirements 
may  be  a limiting  factor.  Secondly,  it  is  difficult  to 
estimate  the  error  introduced  by  imposing  conservation  of 
average  energy  and  stationarity  of  a Maxwellian. 
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Appendix  C 


Numerical  Solution  of  the  Relaxation 
Electron-Electron  Collision  Term 
(Subroutine  Funeer) 


In  this  appendix,  the  numerical  method  used  to 
evaluate  the  relaxation  electron-electron  collision  term 
is  presented.  The  basic  problem  is  defining  a Maxwellian 
such  that  the  term  retains  the  desired  conservation  proper- 
ties. 

The  relaxation  collision  term  is  given  by: 

b>r\  n(S')  - n(Q 


where  n(£)  is  the  Maxwellian  distribution  with  the  same 
average  energy  and  number  density  as  n(£).  Equation  (^2) 
is  easily  cast  into  a finite-difference  form  compatible 
with  the  other  terms  of  the  Boltzmann  equation.  The  result 
is 


b n ; _ n ^ - 0 1 

" ' rce~~ 


(91) 


An  explicit  expression  must  now  be  assumed  for  nL . This 
is  not  necessary  in  evaluating  the  Fokker-Planck  term  where 
only  the  ratio  n^/nj  is  required.  By  analogy  with  the 
continuous  Maxwellian 


where  A and  T are  constants  to  he  determined. 

To  conserve  total  number  density  and  average  energy, 
the  following  requirements  must  be  imposed i 


N Nl 

Z^L^ 

*■->  L=i 


The  sums  on  the  left  of  these  equations  are  recognized  as 
the  total  number  density  and  average  energy  (to  within  a 
factor  of  A £ ) . Therefore 


L-  l 


ne 


:/a£ 


neI/ 


A£ 


(93) 

(94) 


From  Fig.  8 

£;  = (i  — Wat  (95) 

where  A L is  the  energy  interval  used.  Using  Eqs  (92)  and 
(95)*  Eqs  (93)  and  (94)  become 


A«'V'ATI<i-' «YU£/T=  n‘4t 


i=i 

tJ 


A n«E/£ 


L- I 


Ag. 


(96) 

(9?) 


Dividing  Eq  (97)  by  Eq  (96) , the  result  is 


where 


(98) 


“ 3/Z  -L  AE/t 

5 , = i;G-V)  e. 


(99) 


t=  I 


r^-  '/?  ” LAC/-p 


(100) 


To  conserve  average  energy,  the  value  of  T must 
be  such  that  Eq  (98)  is  satisfied.  This  is  accomplished 
using  an  iterative  procedure  with  the  first  estimate  for 
T being  approximately  two- thirds  the  average  energy. 
Subsequently,  the  last  calculated  value  of  T is  used  as  a 
first  estimate.  For  all  calculations  here,  the  value  of 
T is  accepted  when  the  difference  between  the  actual  and 
calculated  average  energy  is  less  than  or  equal  to  10~10. 
Knowing  T,  the  value  of  A is  easily  calculated  from 
either  Eq  (96)  or  Eq  (97).  . . 

The  finite-differenced  Maxwellian  is  now  given  by 
Eq  (92)  and  the  relaxation  term  evaluated  using  Eq  (91). 
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Appendix  D 


Cross  Sections  and  Rates 

The  momentum  transfer  rate  R can  be  expressed  in 
terms  of  the  momentum  transfer  cross  section  as 

MO  = (£)V(o 

where  £ and  m are  the  electron  energy  and  mass  respec- 
tively. Table  I gives  the  momentum  transfer  cross  sections 
used  throughout  this  work.  These  cross  sections  have  been 
constructed  to  yield  a constant  rate.  Therefore,  to  compute 
this  rate  it  is  sufficient  to  consider  only  one  value  of 
the  energy.  For  £ = 1.25  ev 

2 x 1.25  x 1.602  x 10~1Z\l/Z  / 

8.940  x 10 

9.109  x 10_^°  J \ 

and 

Rm  = 5*928  x lO’^cnP-sec"’^  (101) 

Similarly,  the  excitation  rate  R can  be  expressed  as 

R(e)  = (~)  z (o 

Table  II  gives  the  excitation  cross  sections  used  here. 

These  have  been  constructed  to  give  a zero  rate  below 
threshold  and  a constant  rate  Rq  above  threshold.  The 
threshold  energy  for  all  calculations  is  8 ev.  Using  the 


Table  I 


y. 


Momentum-Transfer  Cross  Sections 

8 (ev)  0“m(cm2) 

10.25  3.12x10"16 


£(ev) 

<Tra(om2) 

.25 

2.00xl0“15 

.75 

2.16x10“15 

1.25 

8.94xl0~16 

1.75 

7.56xlO"16 

2.25 

6.67x10"16 

2.75 

6.03xl0_l6 

3-25 

5.55xl0‘16 

3.75 

5.15xlo’16 

4.25 

4.85xlo“16 

4.75 

4.59xl0“16 

5.25 

4.36xl0"16 

5.75 

4.17x10“16 

6.25 

4.00xl0“l6 

6.75 

3,85xlO"16 

7.25 

3.71X10-16 

7.75 

3.59xl0“16 

8.25 

3.48xl0-16 

8.75 

3.38xl0"16 

9.25 

3.29xl0“16 

9.75 

3.20x10"16 

10.75 

3.05xl0'16 

11.25 

2.98xl0-16 

11.75 

2.92xl0-16 

12.25 

2.86xl0_l6 

12.75 

2.80xl0_l6 

13.25 

2.75xl0"16 

13.75 

2.70xl0"16 

14.25 

2.65xlO“16 

14.75 

2.60xl0_l6 

15.25 

2.56xl0"16 

15.75 

2.52xl0"16 

16.25 

2,48xl0"16 

16.75 

2,44xl0"16 

17.25 

2.4lxl0"16 

17.75 

2.37xl0"16 

18.25 

2.34xl0“l6 

18.75 

2.31xl0_l6 

19.25 

2.28xl0‘l6 

19.75 

2,25xl0*16 
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Table  II 

Excitation  Cross  Sections 


£(ev) 

cm  ) 

£(ev) 

crx(cm*) 

.25 

0.00 

10.25 

9.00xl0"17 

.75 

0.00 

10.75 

8.80xl0"17 

1.25 

0.00 

11.25 

8.60xl0”17 

1.75 

0.00 

u.75 

8.40xl0“17 

2.25 

0.00 

12.25 

8.20xl0~17 

2.75 

0.00 

12.75 

8.00xl0"17 

3*25 

0.00 

13.25 

7.90xl0“17 

3-75 

0.00 

13.75 

7.70xl0”17 

4.25 

0.00 

14.25 

7.60x10-17 

4.75 

0.00 

14.75 

7.50x10~17 

5.25 

0.00 

15.25 

7.40x10”17 

5.75 

0.00 

15.75 

7.20xl0~17 

6.25 

0.00 

16.25 

7.10xl0-17 

6.75 

0.00 

16.75 

7.00xl0"17 

7.25 

0.00 

17.25 

6.90X10"17 

7.75 

0.00 

17.75 

6.80X10"17 

8.25 

l.OOxlO"16 

18.25 

6.70X10"17 

8.75 

9.70X10"17 

18.75 

6.60xi0-17 

9.25 

9.40x10’17 

19.25 

6.50X10’17 

9.75 

9.20X10"17 

19.75 

6.50xlo“17 
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Appendix  E 


Analytic  Solution  Using  the  WKB  Method 


Equation  (44)  can  be  solved  using  a combination  of 
the  WKB  method  and  variation  of  parameters.  The  WKB  method 
yields  the  following  solution  to  the  associated  homogeneous 
equation* 


-A"Vu 


4 Cze 


where  and  Cg  are  arbitrary  constants  of  integration. 

Using  variation  of  parameters,  a particular  solution  is 
given  by 


where 


and 


■vfc.fj 


c J- 


a 


Making  the  appropriate  substitutions,  V^y)  and  V2(y) 


may  be  expressed  as 


A">-irc, 


, A'V^V 
s e 

-wrHyy’ 


(103) 


Clearly  V2(y)  = -V^C-y),  and  the  particular  solution  may 


be  expressed  as 


nf (y) = tj  [v, (^  e - v, (- ^ e. 

Therefore,  the  general  solution  to  Eq  (44)  is 

n^  = 'J“L^e  + - V,(-j)e  (104) 

where  the  term  involving  has  been  disregarded.  This 

term  is  unbounded  and  physically  unacceptable. 

Attempts  to  perform  the  integration  in  Eq  (103)  have 
been  unsuccessful.  Furthermore,  if  the  terms  involving 
V^(y)  are  ignored,  the  constant  in  Eq  (104)  cannot  be 
evaluated.  Even  if  this  obstacle  could  be  overcome,  a 
severe  computational  difficulty  is  encountered?  the  constant 
A in  the  exponent  is  generally  a very  large  number. 
Therefore,  the  individual  terms  in  Eq  (104)  are  at  best 
difficult  to  calculate.  For  these  reasons,  a more  straight- 
forward method  is  used  to  solve  Eq  (44)  in  Section  III. 
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Appendix  F | 

Steady-State  Distributions  (Numerical  Results) 

This  appendix  contains  the  numerically  calculated 
steady-state  distribution  functions.  Appendices  B and 
C contain  details  of  the  numerical  procedures.  Plots 
of  the  normalized  distributions,  n(£) , are  presented 
in  Figs.  6 through  8. 


85 


in 

VO 

CO 

o 

OJ 

-it 

VO 

CO 

ON 

Y— 

KA 

in 

0- 

Ov 

Y- 

kv 

in 

K- 

Ov 

Y— 

r— 

*— 

*— 

CVJ 

cvj 

CM 

OJ 

CM 

CVJ 

kv 

KV 

ro 

kv 

ro 

-it 

-it 

^t 

it 

-t 

in 

OJ 

| 

i 

| 

1 

i 

1 

1 

1 

1 

1 

1 

i 

i 

i 

1 

1 

1 

1 

| 

i 

s. 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

KV 

T— 

Y — 

Y — 

Y— 

Y — 

Y— 

Y— 

Y — 

Y— 

Y— 



Y— 

Y 



Y— 

Y— 

Y— 

Y— 

Y— 

1 

X 

X 

X 

w 

X 

X 

X 

X 

X 

X 

M 

X 

X 

X 

X 

X 

X 

X 

X 

X 

t> 

-it 

in 

o 

CVJ 

o 

o 

in 

CTv 

CO 

Y~ 

Cvj 

CO 

co 

-M- 

-it 

CO 

VO 

vo 

VO 

<D 

in 

CVJ 

CO 

in 

UV 

CO 

-t 

-it 

o 

kv 

c*- 

in 

o 

VO 

CM 

CO 

in 

CM 

OV 

la 

CO 

- 

- 

CVJ 

kv 

■3* 

VO 

CO 

- 

- 

CM 

CM 

KV 

kv 

<t 

-it 

in 

VO 

VO 

VO 

Ov 


Q) 

(H 

-O 

a 

Eh 


ir\ 

i 

o 


ii 


CD 

P 

<D 

4-> 

03 

4-> 

CO 

I 

>> 

HD 

a 

<D 

■p 

to 

o 

9 

£ 

■ 

a 

a> 

-X 

o 


II 

Y— 

Y— 

kv 

in 

r- 

Ov 

CM 

-it 

VO 

00 

O 

OJ 

-3-  VI 

1 

OJ 

-it 

vo 

t- 

Ov 

Y— 

Y~ 

y— 

Y— 

y— 

CM 

OJ 

CM 

CM 

OJ 

KV 

fi'N 

KV  ri 

IVJU 

> 

o 

1 

1 

1 

1 

1 

1 

1 

1 

| 

i 

1 

1 

1 

1 

1 

1 

1 

1 

0) 

o 

o 

O 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

O 

•* 

KA 

y— 

Y 

Y 



YT— 

Y 

Y— 

Y— 

Y— 

Y— 

Y— 

Y— 

Y- 

Y— 

Y— 

Y— 

Y— 

Y 

Y— 

CO 

i 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

B 

kv 

Y— 

-it 

Ov 

KV 

CM 

kv 

ro 

CM 

o 

Ov 

o 

-it 

vo 

Y— 

ov 

Ov 

Y 

in 

1 

o 

h- 

Y— 

o 

vo 

OJ 

r- 

• 

Y 

Y 

kv 

VO 

-it 

o 

CM 

in 

K- 

o 

-it 

r- 

o 

s— * 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

a 

CM 

-it 

VO 

CO 

Y— 

CM 

kv 

-it 

in 

vo 

CO 

Y~ 

Y— 

OJ 

OJ 

CVJ 

N. 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

> 

<D 

OJ 

N- 

CM 

h- 

CM 

h- 

CM 

K- 

CM 

K- 

CM 

K- 

CM 

t — 

CM 

CM 

t'- 

CM 

o 

o 

Y- 

Y~ 

OJ* 

CM 

KV 

KV 

-it 

-it 

in 

in 

vo 

vo 

CO 

CO 

Ov 

aJ| 

T— 

Y— 

Y~ 

y — 

Y— 

i-  k 


X—Y. 

Y— 

Y~ 

Y— 

CM 

CM 

CM 

KV 

KV 

KV 

<t 

-it 

-it 

in 

in 

VO 

K- 

Ov 

Y— 

CM 

O 

1 

1 

| 

1 

1 

1 

| 

1 

1 

1 

1 

| 

i 

1 

1 

1 

1 

V 

o 

o 

O 

o 

o 

O 

O 

O 

o 

O 

o 

o 

o 

o 

o 

o 

O 

o 

o 

KV 







Y 

Y— 

Y— 

Y— 

Y 

Y— 

YT— 





Y— 

Y 

Y 

Y— 

Y 



Y 

1 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

t» 

o 

KV 

VO 

ov 

KV 

vo 

in 

o 

OV 

VO 

CM 

CO 

in 

-3* 

KV 

CM 

oo 

<D 

-it 

in 

o 

-it 

VO 

-it 

00 

CM 

-M- 

OV 

OJ 

-it 

-it 

VO 

LT\ 

-it 

vo 

ia 

Y— 

vo 

KV 

- 

vo 

KV 

Y“ 

vo 

KV 

VO* 

KV 

vo 

OJ 

CO 

CM 

KV 

Y— 

KV 

KV 

KV 

ro 

tn 

CM 

CM 

0J 

Y- 

Y- 

o 

O 

o 

| 



Y— 

Y— 

Y 

Y~ 

Y — 

Y— 

Y— 

Y 

Y— 

Y— 

Y— 

Y 

ov 

o 

K- 

in 

-it 

> 

o 

o 

O 

O 

o 

O 

o 

O 

o 

O 

o 

O 

O 

O 

O 

o 

O 

o 

o 

CD 



Y— 

Y— 



Y— • 

Y 

Y 

r— 

Y— 

Y— 

Y 

Y— 

Y— 

Y— 

Y 

Y 

Y— 

Y — 

T— 

KV 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

1 

CO 

vo 

CM 

Ov 

o 

Ov 

-it 

KV 

o 

vo 

OV 

CM 

N- 

CO 

ro 

TO 

o 

vo 

OJ 

0 

Ov 

vo 

-it 

CO 

o 

vo 

KV 

VO 

OJ 

in 

K- 

VO 

VO 

KV 

00 

Y- 

a 

VO 

LTV 

KV 

-- 

uv 

CM 

vo 

KV 

*- 

KV 

- 

Is- 

CM 

-it 

vd 

Y— 

> 

0) 

u> 


OJ  K- 


in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

uv 

in 

CM 

K- 

CM 

K- 

OJ 

h- 

OJ 

K- 

CM 

K- 

CM 

K- 

CM 

K- 

OJ 

h- 

CM 

c- 

CM 

CM 

KV 

KV 

-it 

-it 

in 

in 

vo 

vo 

h- 

CO 

CO 

Ov 

ov 

86 


77x10  5.63x10“  J 19.75  3.09xiO"JU  6.96 


t>  in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

m 

in 

in 

m 

in 

in 

in 

<U  CM 

n- 

CM 

r- 

CM 

t'- 

CM 

n- 

CM 

c- 

CM 

c- 

CM 

c- 

CM 

f- 

CM 

t-- 

CM 

<J 0 

CM 

CM 

m 

m 

*3” 

Ln 

in 

VO 

VO 

N- 

co 

CO 

CT\ 

8? 


9.75 


Table  V 

-Planck  Steady-State  Dist.  -n«/n  = lO-5,^  = 10"18,  £ = .732 


o 

Y- 

Ol 

OJ 

in 

NO 

NO 

K- 

CO 

ON 

On 

O 

Y- 

Y- 

CM 

KN 

<J" 

▼— 

T— 

V— 

y— 

T“ 

T— 

▼— 

Y“ 

r— 

V— 

Y— 

Y— 

o — 

Y— 

CNJ 

Ol 

OJ 

OJ 

OJ 

1 

CM 

OJ 

i 

1 

1 

i 

1 

1 

1 

i 

i 

1 

1 

1 

1 

1 

1 

1 

1 

1 

V 

o 

O 

o 

o 

o 

O 

o 

o 

o 

o 

o 

O 

o 

O 

O 

o 

o 

o 

o 

O 

Y — 

Y— 

Y— 

Y— 

T— 

T— 

Y— 

Y— 

y — 

Y— 

Y— 

1 — 

Y— 

T— 

Y— 

Y— 

Y— 

Y— 

Y 

r— 

| 

* 

X 

X 

X 

X 

X 

X 

X 

x 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

> 

in 

VO 

ON 

VO 

K- 

ON 

ON 

VO 

o 

Y — 

K> 

On 

o 

ro 

O') 

CO 

CO 

o 

rO 

ON 

(D 

OJ 

VO 

in 

ON 

ON 

ON 

UN 

OJ 

Ol 

kn 

VO 

NO 

O 

c- 

o 

o 

t- 

la 

OJ 

On 

Y~* 

KN 

t— 

J- 

KN 

NO 

Y— 

Ol 

*3* 

CO 

Y— 

in 

tn 

OJ 

in 

t^- 

i 

t> 

NO 

NO 

UN 

-3" 

in 

Ol 

CM 

o 

1 

1 

OJ 

in 

i 

in 

i 

1 

UN 

1 

VO 

I 

VO 

1 

t- 

1 

<D 

o 

o 

O 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

m 

i — 

T— 

Y— 

Y— 

T— 

Y— 

Y — 

Y— 

Y — 

Y — 

Y— 

Y— 

Y— 

Y— 

Y— 

Y— 

v— 

Y 

Y— 

i 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

B 

On 

K's 

Ol 

Y— 

On 

UN 

ON 

l- 

-=r 

-3* 

T— 

CO 

X- 

CO 

CO 

rO 

o 

NO 

ro 

KN 

O 

UN 

CM 

f- 

ro 

00 

K- 

Y— 

VO 

Y— 

c- 

VO 

Ol 

->T 

NO 

NO 

NO 

in 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

q 

K- 

KN 

NO 

OJ 

in 

C\l 

-3- 

ON 

Y~ 

NO 

Y— 

CM 

CO 

Y- 

KN 

. — > 

UN 

UN 

UN 

UN 

UN 

UN 

in 

UN 

in 

UN 

UN 

in 

U 

UN 

in 

m 

UN 

UN 

UN 

i> 

OJ 

K- 

Ol 

K- 

Ol 

K- 

CM 

OJ 

K- 

CM 

CM 

K- 

CV1 

C— 

OJ 

t— 

CM 

O 

O 

Y- 

CM 

OJ 

KN 

<r 

<1- 

UN 

UN 

NO 

NO 

t^- 

00 

CO 

ON 

uJ 

T~ 

Y— 

Y— 

Y— 

Y— 

Y— 

Y~" 

Y— 

Y— 

Y“ 

Y— 

Y— 

Y— 

Y— 

Y“ 

Y~ 

Y— 

Y- 

Y- 

CM 

CM 

OJ 

KN 

UN 

m 

UN 

NO 

NO 

X- 

f— 

00 

ON 

On 

OJ 

o 

1 

1 

1 

1 

1 

1 

i 

\ 

1 

1 

■ 

1 

1 

1 

1 

1 

1 

1 

r 

o 

o 

o 

o 

O 

o 

O 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 



1 — 

Y— 

Y 



Y 

Y — * 

Y— 

Y 

Y 

Y— 

Y 

Y 

Y 

Y 

Y— 

Y— 

Y— 

Y— 

Y 

1 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

> 

-3- 

o 

CM 

ro 

u 

NO 

-it 

UN 

ON 

o 

ON 

ON 

Y— 

o 

in 

o 

ON 

NO 

NO 

Is- 

CD 

CO 

CO 

m 

KN 

CM 

K- 

h- 

*3- 

00 

OJ 

CM 

UN 

o 

o 

in 

o 

o 

- 

NO 

CM 

ON 

ro 

Y— 

NO 

CM 

CO 

— 

-=t 

Y— 

UN 

CM 

UN 

-x 

o 


UN 

in 

in 

UN 

-=t 

KN 

KN 

KN 

CM 

OJ 

^ 

T-  O 

O 

1 



Y — 

Y— 

Y— 

Y 



Y— 

Y 

Y— 

Y— 

Y— 

Y 

Y 

Y—  Y— 

Y— 

On 

co 

CO 

o 

o 

o 

o 

O 

O 

O 

O 

O 

O 

o 

O 

O 

O O 

o 

O 

o 

o 

o 

o 

Y 

Y 

Y— 

Y 



Y 

Y 

Y— 

Y — ’ 

Y— 

Y 

T— 

Y—  Y— ' 

Y— 

Y— 

Y— 

Y— 

T— 

KN 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X X 

X 

X 

X 

X 

1 

On 

On 

CM 

ro 

CO 

CM 

"3- 

o 

-^r 

co 

CM 

-J-  CO 

ON 

CO 

UN 

n 

1 

Y~ 

C0 

00 

OJ 

Y— 

T“ 

UN 

CM 

o 

CO 

O 

Y-  Y- 

r^\ 

ON 

UN 

KN 

On 

UN 

CM 

- 

UN 

OJ 

CO 

K> 

UN 

CM 

r<> 

— <)■ 

- 

in 

vo 

t> 

UN 

in 

UN 

in 

m 

UN 

UN 

in 

Lf\ 

UN 

in 

UN 

in 

UN 

in 

in 

in 

UN 

in 

(D 

Ol 

K- 

CM 

K- 

CM 

CM 

t- 

OJ 

K- 

OJ 

K- 

CM 

K- 

CM 

K- 

CM 

t- 

CM 

oO 

- 

- 

OJ 

CM 

UN 

in 

VO 

vo 

co 

CO 

On 

88 


19.75 


.D-A064  212 


INCLASSIF1ED 


AIR  FORCE  INST  OF  TECH  MRIGHT-PATTERSON  AFB  OHIO  SCH — ETC  F/G  20/8 
ELECTRON-ELECTRON  COLLISIONS,  a SIMPLE  MODEL  FOP  PARTIALLY  IONI— ETC(U) 
DEC  78  ML  HASTINGS 


1 


AFIT/6EP/PH/78D-6 


NL 


in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

m 

in 

in 

in 

in 

in 

in 

in 

in 

in 

CVJ 

h- 

CM 

t*- 

CM 

CM 

h- 

CM 

t'- 

CM 

h- 

CM 

h- 

CM 

t*- 

CM 

r- 

CM 

t'- 

*- 

•r- 

cvi 

CM 

in 

in 

-sr 

in 

in 

VO 

ID 

oo 

OO 

ON 

ON 

-sc-znxrr! r 


cu 

s. 

K\ 

I 

t> 

«>i 


in 

in 

in 

in 

in 

in 

VO 

VO 

VO 

VO 

VO 

c- 

t- 

h- 

h- 

OO  CO  CO 

1 

1 

1 

1 

i 

i 

1 

1 

i 

i 

1 

1 

1 

i 

i 

1 

1 

1 

1 

| 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

T— 

▼— 

•— 

*— 

▼— 

v— 

V— 

T— 

, — 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

x 

x 

X 

h- 

CVJ 

h- 

m 

CO 

00 

r— 

o 

ov 

CO 

o 

in 

Ov 

00 

CVJ 

CVJ 

h- 

co 

CO 

CM 

m 

in 

o 

f- 

in 

o 

Ov 

tn 

oo 

00 

CO 

in 

t-- 

VO 

N- 

b- 

in 

m 

CM 

«- 

VO 

-Nf 

in 

oo 

in 

m 

CM 

t'- 

in 

• 

in 

o 

oo 


It 

IUI 

9k 

VO 

T— 

I 

0 
II 

A 

in 

1 

o 


M || 
M o 
M .C 

t>  > 

C 

<U  I 

H 

& . 
a *> 
eh  n 


I 

s> 

(U 

tn 

i 

B 

o 


«-  u) 


CVJCVJCVJ*—  *r-000000 

oooooooooooo  o ’"o  *"o  °o  °o  °o  °o  °o 

XXXXXXXXXXXXXXXXXXX^ 

cvj^-vo^-N-mt^-^-fnfOvinoNvoovb-mcnob-co 

minocucoincvjvo 

w <-  >-  n in  tn  ai  ^ r-  f-  m ^ cl  r-‘  t-  </  in  cJ 


in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

Lf\ 

m 

in 

CVJ 

C'- 

CVJ 

CVJ 

t'- 

CVJ 

t- 

CVJ 

c- 

OJ 

o- 

CVJ 

t- 

CVJ 

t*- 

cvj 

c- 

CVJ 

t- 

o 

o* 

OJ 

cvj 

tn 

in 

*3* 

T— 

-s- 

in 

in 

VO 

VO 

t*- 

t- 

co 

co 

• 

Ov 

• 

Ov 

<u 

4» 

s 

CO 

I 

N 

•a 

as 

<u 

+> 

CO 

§ 

•H 

p 

g 

«fl 

« 


OJ 

K\ 

I 

> 
a> 
' — - 

in 


x 

o 


i i i i i i 

o o o o o o 

X X X X >5  X 


VO  <f 


CVJ 

CVJ 

CVJ 

CVJ 

m 

in 

tn 

tn 

m 

•3* 

1 

i 

1 

1 

i 

i 

i 

i 

i 

i 

i 

1 

1 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

*■“ 

T— 

*■“ 

T— 

T— 

T— 

*— 

T— 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

00 

CO 

CVJ 

in 

in 

o 

OV 

tn 

t- 

o 

f- 

CVJ 

ov 

oo 

in 

t"- 

vo 

tn 

cvj 

tn 

CVJ 

CO 

CO 

tn 

cvj 

T— 

- 

in 

tn 

OJ 

in 

CVJ 

*>! 

0) 

tn 

i 

Bj 

CJ 


o o 

x x 
in 


n 

in 

in 

in 

in 

•St* 

<1- 

<r 

tn 

tn 

tn 

tn 

CVJ 

CVJ 

CVJ 

*■” 

*■" 

t— 

T— 

r— 

▼— 

T— 

*— 

T— 

o 

o 

o 

o 

o o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

T™ 

*■" 

T— 

T— 

v- 

T— 

r— 

X 

X 

X 

X 

X X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

o 

o 

vo 

in 

in  — 

tn 

in 

tn 

h- 

y— 

o 

tn 

Ov 

OV 

ro 

OJ 

vo 

ov 

o in 

rn 

o 

vo 

CO 

CVJ 

00 

o 

co 

<1- 

•3* 

in 

t*- 

• 

• • 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

m 

CVJ 

— h- 

in 

m 

CVJ 

oo 

vo 

tn 

oo 

in 

tn 

> 

in 

in 

in 

in 

in 

in 

in 

in 

in 

m 

in 

m 

in 

in 

in 

in 

in 

in 

in 

in 

4) 

CVJ 

• 

• 

cvj 

• 

• 

CVJ 

• 

• 

CVJ 

• 

h- 

• 

CVJ 

• 

• 

CVJ 

• 

• 

CVJ 

• 

• 

cvj 

• 

i^- 

• 

CVJ 

• 

• 

OJ 

• 

h- 

• 

at 

CVJ 

CVJ 

tn 

tn 

in 

in 

vo 

vo 

t- 

00 

oo 

Ov 

Ov 

91 


_ 


Harold  L.  Hastings  was  born  on  16  April  19^9  in 
Dallas,  Texas.  He  graduated  from  W.  W.  Samuell  High 
School  in  Dallas,  Texas,  in  1967.  In  June  1971,  he 
received  the  degree  of  Bachelor  of  Science  in  physics 
from  Stevens  Institute  of  Technology,  Hoboken,  New 
Jersey.  Upon  graduation,  he  entered  the  Officer 
Training  School  and  received  a commission  as  a Second 
Lieutenant  in  the  USAF.  He  then  entered  Undergraduate 
Pilot  Training  at  Webb  APB,  Big  Spring,  Texas.  He 
served  as  co-pilot  and  aircraft  commander  in  C- 130's 
until  entering  the  School  of  Qigineering,  Air  Force 
Institute  of  Technology,  in  June  1977. 


Permanent  address:  60  E.  Westfield  Avenue 


Roselle  Park,  New  Jersey 


1 


UNCLASSIFIED, 


REPORT  DOCUMENTATION  PAGE 

^ ’ - 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER  2.  GOVT  ACCESSION  NO. 

AFIT/GEP/PH/78D-6 

3.  RECIPIENT'S  CATALOG  DUMBER 

4.  TITLE  (and  Subtitle) 

ELECTRON-ELECTRON  COLLISIONS! 

A SIMPLE  MODEL  FOR 

PARTIALLY  IONIZED  GASES 

5.  TYPE  OF  REPORT  A PERIOD  COVERED 

MS  THESIS 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHORfa; 

Harold  L.  Hastings 

Capt  USAF 

- 

8.  CONTRACT  OR  GRANT  NUMBERfsJ 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS  ^ 

Air  Force  Institute  of  Technology 
(AFIT-EN)  , , 

Wright-Patterson  AFB,  Ohio  45433 

10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  4 WORK  UNIT  NUMBERS 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Air  Force  Weapons  Laboratory 
(AFWL/DYP) 

K^rtlanri  AFRr  New  Mexico  87117 

12.  REPORT  DATE 

December,  1978 

13.  NUMBER  OF  PAGES 

100 

14.  MONITORING  AGENCY  NAME  & ADDRESS^//  dilferent  Irom  Controlling  Olticc) 

15.  SECURITY  CLASS,  (of  this  report) 

UNCLASSIFIED 

15a.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (ol  this  Report) 

Approved  for  public  release i distribution  unlimited 

17.  DISTRIBUTION  STATEMENT  (ol  the  abstract  entered  in  Block  20.  II  dillerent  Itom  Report) 

18.  SUPPLEMENTARY  NOTES 

Approved  for  public  release}  IAW  AFR  190-17 

J.  pA?pp'^Tm^,  USAF 

Director  of  Information  /"<*  > '// 

19.  KEY  WORDS  (Continue  on  reverse  side  if  necessary  and  identify  by  block  number) 

Electron-Electron  Collisions 

Fokker-Planck  Equation 

Partially  Ionized  Gases 

\\  

20.  ATTRACT  (Continue  on  reverse  side  If  necessary  and  Identify  by  block  number)  . , . 

The  electron-electron  collision  term  of  the  steady-state 
Boltzmann  equation  is  replaced  by  a simple  relaxation  term  that 
effectively  linearizes  the  equation.  Additional  assumptions  are 
made  to  simplify  the  equation  further.  Ionization  is  ignored  and 
a two-level  atom  is  assumed.  The  effects  of  a DC  electric  field 
are  included.  Using  the  relaxation  term  to  account  for  electron- 
electron  collisions,  an  approximate  analytic  solution  is  derived. 
Temporal  and  steady-state  characteristics  of  the  relaxation  termer 

DD 


FORM 
AN  73 


1473 


EDITION  OF  1 NOV  65  IS  OBSOLETE 


_UNCLA.SSIFJ.ED_ 


/’’V 

l 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (H?icn  Dm*  Entered) 


CURITY  CLASSIFICATION  OF  THIS  P AGE  (H?»an  DmU » Entmrmd) 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  Tu.«“  P A5Cf  M?i*n  Pmtm  kntermd) 


